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ABSTRACT 


Subsonic  Kernel  function,  transonic  box,  and  supersonic  box  methods 
for  computing  unsteady  aerodynamics  are  applied  to  the  problem  of  inter¬ 
action  of  a  general  trapezoidal  wing  with  a  downstream  rectangular  control 
surface  lying  in  the  wake  of  the  wing.  The  unsteady  aerodynamic  forces 
are  related  to  a  set  of  collocation  stations  through  a  series  of  matrix 
transformations,  interpolations,  and  differentiations.  The  resulting 
matrix  is  a  set  of  aerodynamic  influence  coefficients  (AICs)  that  are 
directly  applicable  to  flutter  analysis. 

The  transformation  of  the  unsteady  aerodynamics  into  AICs  is  pre¬ 
sented  as  a  separate  discussion;  followed  by  discussions  for  the  develop¬ 
ments  of  analytical  techniques  for  each  flight  regime.  The  analytical 
developments  and  a  discussion  of  the  basic  single-planar-surface  are 
presented,  followed  by  the  complete  two-surface  solutions  for  the  general 
aerodynamic  forces.  Each  of  the  three  numerical  methods  is  developed  by 
detailing  the  complete  set  of  equations  necessary  to  compute  airloads  on 
the  configurations  considered.  A  computer  program  to  determine  the  AIC 
matrix  for  each  flight  regime  is  presented  with  a  complete  discussion  of 
usage  and  logical  flow.  Also  included  are  program  listings,  flow  charts 
and  sample  input,  and  output  problems. 


PART  I 


INTRODUCTION 


The  requirements  for  determination  of  flutter  margins  of  safety 
for  the  lifting  surfaces  of  .  advanced •  guided  missiles  have  precipitated 
a  need  for  accurate  methods  of  analysis  of  unsteady  aerodynamic  loading 
in  the  high  subsonic,  transonic,  and  supersonic  flight  regimes.  These 
methods  must  not  only  account  for  the  high  degree  of  chordwise  and  span- 
wise  deformation  of  the  surfaces,  but  also  include  the  interference 
effects  between  tandem  lifting  surfaces.  Recent  developments  in  lifting 
surface  theory  in  the  three  Mach  number  regimes  have  permitted  extensions 
(Refs.  2  and  3)  to  determine  airloads  on  typical  missile  wings  with  down¬ 
stream  control  surfaces.  These  extensions  account  for  the  interaction 
and  wake  effects  as  well  as  for  the  three-dimensionality  of  the  flow 
for  a  trapezoidal  wing  planform  and  a  coplanar  rectangular  control 
surface  placed  at  an  arbitrary  distance  downstream  of  the  unswept  trailing 
edge  of  the  wing.  An  underlying  assumption  in  these  methods  is  that  the 
missile  body  diameter  is  large  enough  compared  to  the  spans  of  the  surfaces 
that  the  body  surface  acts  as  a  reflection  plane  for  distrubances  at  the 
line  of  its  iirtersection  with  the  lifting  surfaces. 

The  present  study  extends  the  methods'  of  Ref.  3  to  obtain  aerodynamic 
influence  coefficients  (AIC's)  that  relate  the  forces  on  the  surfaces  at 
a  discrete  set  of  points  (control  or  collocation  points)  to  the  tranverse 
deflections  of  the  same  set  of  points.  Subsonic  kernel  function,  transonic 
box,  and  supersonic  box  methods  for  computing  the  oscillatory  AIC's  are 
applied  to  the  interference  problem  of  a  general  trapezoidal  wing  with  a 
downstream  rectangular  control  surface  lying  in  the  wake  of  the  wing. 

Highly  efficient  numerical  methods  for  computation  by  the  kernel  function 
and  Mneh  box  techniques  have  been  employed,  along  with  the  techniques  for 
the  newly  developed  transonic  box  method,  to  obtain  AIC's  which  account 
for  all  interference  effects  within  the  framework  of  linearized  theory. 


Discussion  of  the  derivation  of  the  AIC's  is  given  in  Part  III  for 
the  three  Mach  number  regimes.  The  analytical  basis  of  the  theories  are 
outlined  in  Appendices  to  Part  III.  Each  of  the  three  numerical  methods 
is  discussed  and  the  basic  equations  necessary  to  compute  airloads  on 
the  configurations  considered  are  summarized. 

The  three  computer  programs  are  presented  in  Parts  IV,  V,  and  VI, 
for  the  subsonic,  sonic,  and  supersonic  cases,  respectively.  In  addition 
to  a  technical  outline  of  the  methods  employed,  each  Part  is  a  manual 
containing  a  complete  discussion  of  usage  and  logical  flow  accompanied 
by  program  listings,  flow  charts  and  sample  input  and  output.  Each  Part 
also  presents  results  confuted  by  operation  of  the  program  for  typical 
planforms. 
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PART  II  -  NOMENCLATURE 


< 

Free-stream  acoustic  velocity 

Coefficients  of  Kernel  function  pressure  series 

L  cal  semi -chord 

Reference  semi-chord 

Element  of  dimensionless  AIC  matrix 

Pressure  coefficient 

Element  of  Kernel  function  matrix 


Supersonic  source  influence  function 
AIC  control  point  displacement. 


Kernel  function 
Local  reduced  frequency 


Reference  reduced  frequency 


Mach  number 


Complex  amplitude  of  pressure 
Planform  area 


Semi-span 

Free-stream  fluid  velocity 
Chebyshev  polynomial 
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Substantial  derivative  matrix 


Complex  amplitude  of  dovnwash 


X,y,z  Cartesian  coordinate  system 

p2  1-J^  for  M<  lj  J^-l  for  H>1 

A  Box  length 

Cartesian  coordinate  system  variables 

f  .  Special  coordinates  for  collocation  and  integration  points 


Complex  amplitude  of  velocity  potential 
Transonic  doublet  influence  function 


Atmospheric  density 

Angular  frequency  of  harmonic  oscillation  ~  rad/sec 
2 

Aspect  ratio  ~  2s  /S 
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PAST  III 

DISCUSSION  OF  THE  DERIVATION  OF  AERODYNAMIC 
INFLUENCE  COEFFICIENTS 


DERIVATION  OF  SUBSONIC  AERODYNAMIC  INFLUENCE 
COEFFICIENTS 

The  subsonic  kernel  function  procedure  was  developed  by  Watkins, 

Woolston  and  Cunningham/  extended  to  a  wing-tail  combination  by 
2  3 

Moore  and  Park,  and  refined  by  Andrew.  The  present  study  extends 
the  method  of  Ref.  3  to  obtain  aerodynamic  influence  coefficients  (AICs). 
Oscillatory  AICs  have  been  defined  in  Ref.  4  to  relate  control  point 
forces  to  control  point  deflections  through  the  matrix  equation. 

{F)  .  pu2br29[Ch]<h}  (3.1) 

The  derivation  of  the  AICs  requires  a  review  of  the  technique  of  Refs. 

1-3.  (NOTE:  These  references  are  outlined  in  Part  IV,  Section  A).  The 

starting  point  is  an  assumed  series  for  the  lifting  pressure  distribution. 

3 

It  is  chosen  in  the  form 

\!~  2  ~~2  ^  M 

V^>  *  1  I  pm<a>  (3-2) 

n=0  m=0 

i 

where  flie  rhordwise  pressure  functions  arc 

f0(£)  =  \Jv  -|)/(1  +  £>  (3.3a) 

<„(£)  =  V1  -  e 2  un.i<e>'  1 *  «  <3.3b> 

the  spanwise  pressure  functions  are 


(3.4a) 
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(3.4b) 


pmia>  *  a2  om.i«a>.  i  *  "> 

and  the  Chebyshev  polynomial  recurrence  relation  is 

Uk(x)  =  2xUk  l{x)  -  Uk  2(x)  (3.5) 


when  Uq(x)  =  1.0  and  Uj(x)  =  2x.  The  numerical  procedures  of 
Refs.  1  -  3  lead  to  a  matrix  formulation  of  the  aerodynamic  lifting  sur¬ 
face  integral  equation  as  an  equation  between  the  control  point  down- 
washes  and  the  amplitudes  of  the  assumed  pressure  modes 

{w/n}  =  [  D  ]  {a  }  n 

1  w  1  nmJ  1  nm 


The  solution  of  £q.  3.6  by  least-squares  methods  is  given  by  Ref.  3 
which  we  write  in  the  form 

Km1  "  [AM<!]fw/U)  (3-7) 

In  order  to  find  the  AICs  it  is  necessary  to  define  a  grid  of  control 
points.  We  find  it  convenient  to  divide  the  surface  into  NS  strips  and  to 
locate  the  control  points  on  the  centerline  of  each  strip.  We  further 
choose  NC  control  points  on  each  strip  located  at  the  same  fractional 
chord  location  on  each  strip.  A  typical  distribution  of  control  points  is 
shown  in  Figure  3.1  for  a  wing-tail  combination. 

The  downwashes  required  in  Eq.  3.7  are  obtained  by  a  substantial 
differentiation  of  the  AIC  control  point  deflections 

{w/u}  =  f  W]  {h}  (3,8) 

Sim  e  the  downwashes  required  in  Ref.  3  are  at  an  optimum  set  of  points 
different  from  the  AIC  control  points,  the  matrix  [W]  is  seen  to  be  an 
interpolation  as  well  as  a  substantial  differentiation  matrix. 
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DISTRIBUTION  OF  CONTROL  POINTS 


CONTROL 

WING  SURFACE 


Figure  3.1 
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Ftos  the  pressure  raodes  end  their  sag&iiudeg  the  forces  say  be 
found  by  an  integration  procedure*  The  3*2  forces  on  each  strip  nay 
be  found  by  integration  of  the  pressure  on  the  strip  in  the  region 
of  each  of  the  HC  forces  to  obtain  equivalent  concentrated  forces. 
This  leads  to  a  relationship  between  the  forces  and  the  pressure  node 
amplitudes 


{F)  =  <i«2(8]{anm} 


(3.9) 


We  designate  the  matrix  (B]  as  the  integration  matrix. 

Combining  Eqs.  3*7,  3,3  and 3.9,  and  identifying  the  result  with 
Eq.  3.1  leads  to  the  subsonic  oscillatory  AlCs. 

(Ch)  =  (»/2k*)[B)[AM<i][W]  (3.10) 


Ref.  3  provides  j];  the  extension  to  obtain  AICs  requires  the  de¬ 

velopment  of  the  matrices  [W]  and  [B],  These  are  discussed  next. 

THE  SUBSTANTIAL  DIFFERENTIATION  MATRIX  [w] 

The  substantial  differentiation  matrix  is  derived  by  surface  fitting 
techniques.  For  maximum  accuracy  we  fit  the  surface  "in-the-small," 
i.  e.  ,  locally.  Rather  than  use  a  surface  fit  per  se,  we  shall  fit  curves 
in  the  chordwise  direction  and  then  fit  similar  curves  spanwise  along 
lines  of  constant  chord  fraction.  A  higher  order  polynomial  is  not  well 
behaved  between  points,  so  we  choose  to  connect  a  series  of  parabolas. 
A  number  of  options  exists  as  the  number  of  points  is  increased  so  it  is 
well  to  develop  the  equations  systematically. 

With  two  control  points  the  curve,  of  course,  is  a  straight  line, 

Um  equation  may  be  written  in  matrix  form. 


h(x) 


U  xj 


(3.11) 
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With  three  control  points  the  parabola  is  given  by 


(3.12) 


In  the  case  of  four  points  we  fit  two  parabolas  that  are  tangent 
half-way  between  the  second  and  third  points.  For  x  <  x^  ^  -  (xz  +  x^)/2 
we  write 


h(x) 


ao  + 


a  jX 


a2x 


(3.13) 


and  for  x  ->•  x 


2-3 


h(x)  =  bg  +  bjX  4  b^x 


(3.14) 


At  x  ;  j  the  deflections  and  slopes  must  be  equal. 


a0  '  alX2-3  f  a2x2-3  =  b0  +  blx2-3  +  b2x2-3 


(3.15) 


^  ^a2  x2-3  ~  bl  ^  ^b2x2-3  (3.16) 

Tlu;  unknown  coefficients  are  determined  by  the  solution  of  the  matrix 
equation 
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(3.17) 


Solving  for  the  coefficients  leads  to  the  equations  for  the  curve. 


h(x) 


ra  ^ 
a0 


[l  x  x2  0  0  oj( 


>  for 


U1 

b2 

V  J 

''aj 
0 


2-3 


[o  0  0  1  X  X2J( 


b. 


/  for  x  > 


'i 


b2 

V.  J 


l2-3 


(3.18a) 


(3.18b) 


The  cii.Hc  of  five  points  leads  to  the  general  pattern.  We  use  two 
p.t  r.i  l»o l.ts  at  the  ends  as  in  the  four  point  case  and  one  intermediate 
parabola  th.tt  goes  through  one  point  and  is  tangent  to  the  other  parabolas. 
I- or  x  '  (x^  i  x.j)/2 
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h(x)  =  aQ  +  ajx  +  a2x 

(3.19) 

h  '(x)  =  a  j’ '+  Z&2  x 

(3.20) 

For  x2_,{  ^  x  <-  x3_4  -  (x3  i  x4)/2 

h(x)  =  bQ  +  b1  x  +  b2  x2 

(3.21) 

h  '{x)  =  b  ^  +  2b2  x 

(3.22) 

Finally,  for  x  >  x3  4 

2 

h(x)  =  Cq  +  CjX  +  c2x 

(3.23) 

h  '(x)  =  c  j  +  2c2  x 

(3.24) 

The  nine  coefficients  are  found  from  the  solution  of 


The  general  pattern  is  indicated  by  the  partitioning.  The  equations  for 
six  or  more  points  may  be  written  by  inspection.  The  equation  for  the 
curve  may  be  written  in  the  form 


a2 


. 

r  2 

h(x  *-  *2_3)  j 

Ixx  OOO  000 

bl 

h*x2-3  ^  x  ^  *3.4)  1  = 

0  0  0  1  x  x2  0  0  0 

< 

b2 

h(x  2r  x3  4)  J 

000  000  Oxx2 

- 

b3 

C1 


(3,26) 


The  curves  that  are  obtained  from  this  procedure  resemble  the 
deflection  curve  of  a  beam  over  multiple  deflected  supports.  If  only  one 
support  is  deflected  the  curve  damps  out  as  the  distance  from  the  de¬ 
flected  support  is  increased.  Some  examples  of  the  surface  fits  are 
shown  in  Figures  -3,2  -  3.4. 

The  foregoing  procedure  may  be  summarized  formally  by  re¬ 
writing  Kq.  3.2b  as 


f  T(xa)]  {an)  -  [B(xa)]  (ha) 


(3.27) 


and  Kq.  3.26  as 


{h(x)}  ~  [ C(x)l  (an)  (3.28) 

The  deflection  at  an  arbitrary  location  is  found  in  terms  of  the  AIC 
control  point  deflections  by  combining  Eqs.3.27  and  3.28. 
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Figure  3.2 
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TYPICAL  (W)  MATRIX  DEFLECTION  PATTERN 


FOR  A  UNIT  DISPLACEMENT  AT  B 


TYPICAL  (U)  MATRIX  DEFLECTION  PATTERN 


FOR  A  UNIT  DISPLACEMENT  AT  C 


(3.29) 


{h(x)}  =  [C{x)3  [T(xa)3-1  [B(xa)3{ha} 

The  flown  wash  is  the  substantial  derivative  of  Eq.  3.29. 

r  Ifc  '  1  u  (3.30a) 

=  (h{x)}  +  i  ^  {h(x)}  (3.30b) 

r 

The  derivative  of  Eq.3.29  is 

£  {M*)3  =  fC'(x)3  [T(xa)3_1  [B(xa)3{ha3  (3.31) 

where  [C'(x)3  is  the  matrix  of  derivatives  of  C(x). 

We  next  consider  the  interpolation  in  the  spanwise  direction.  We 
have  already  defined  {h  }  as  the  set  of  AIC  control  point  deflections. 

Since  the  collocation  points  in  the  lifting  surface  theory  have  different 
chordwise  and  spanwise  locations,  we  introduce  the  matrix  {h^ }  of  de¬ 
flections  on  the  AIC  strip  centerline  at  the  fractional  chord  locations  of 
the  lifting  surface  theory  collocation  points.  Interpolation  among  the 
{h^}  on  all  of  the  strips  then  leads  to  the  lifting  surface  theory  colloca¬ 
tion  point  deflections  { hc } .  The  substantial  derivative  matrix  is  then 
defined  by 

,W^{ha}  =  S  {hc}  +  1  IT  *hc}  0.32) 

r 

The  chordwise  interpolation  leads  to 

{h()}  fC(xb)l  [T(xa)T1  fB(x;i)l  { h;i }  (3.33) 

•i  ml 

•  |C'(xb)]  I  T(xa)3  “ 1  |B(xa)l  {hn)  (3.34) 
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Then  the  spanwise  interpolation  leads  to 


{hc>  =  [C(yc)]  [T(yair1  [B(ya>]  {hbJ 

(3.35) 

and 

=  [Cfy^HTlyjr'rBly^]^  {hb) 

(3.36) 

The  required  matrix  fW]  =  +  i  f  Wj]  is  found  by  combining  Eqs. 

3.32  -3.36.  The  real  and  imaginary  parts  appear  formally  as 

'V 

=  [C(yc)]  ]T(ya)]_1  [B{ya)l  [C*^)]  [Tfx^]"1  fB(xa)] 

(3.37) 

and 

[Wjj  = 

(kr/br)  [C{yc)]  [  T(ya)]  "1  [B(ya)]  [C{xb]  [Tfxj]"1  [Bfxj] 

(3.38) 

The  formalism  may  be  illustrated  simply  by  considering  two  strips  with 
4  AIC  points  and  4  aerodynamic  collocation  points  as  shown  in  Fig¬ 
ure  3.5. 


Figure  3.5 
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Combining  the  deflections  yields  a  partitioned  form 


The  spanwise  interpolation  on  the  two  spanwise  lines  gives 


on  the  forward  line  and 


h  ■ 
c2 
h 

'  M 

fh, 

b2 

h 

c4, 

b4 

on  the  aft  line.  Combining  the  deflections  leads  to  another  partitioned 
form  but  also  requires  a  row  rearrangement  matrix  to  order  the  {h^J 
properly 


r 


i 


fh 

_ 

”1  0  0  0“ 

rK 

C1 

1*13  !  0 

bi 

h 

1  ~>  1 

0  0  10 

h. 

C2 

1 

b2 

h 

0  10  0 

hb 

c3 

0  !  Yi3 

b3 

hc 

L  !  1 41 

0  0  0  1 

hb 

4. 

4J 

The  same  procedure  is  used  for  the  slopes.  We  note  that  Eqs,  3.37  and 
3.38  must  be  generalized  to  include  the  row  rearrangement  matrix,  and 
the  formal  equations  for  the  real  and  imaginary  parts  of  the  substantial 
differentiation  matrix  become 

N  *  K>]  K»>1  ■'  K»>1  M  Ic,<^>l  hvl  '* 


and 

|Wj|  .  (kr/br)  jc(yc)j  [Tl^))-1  |B(ya)|  |r)  |c(^)|  |T(*a)j-‘  |B(xa)| 

whore  the  format  of  each  factor  is  illustrated  in  the  above  example. 

THE  INTEGRATION  MATRIX  [B] 

The  integration  matrix  converts  the  pressure  distribution  on 
«*neh  ntrlp  into  nn  equivalent  system  of  concentrated  forces  at  the 
A1C  control  points  by  integrating  the  pressure  spsnwise  and  chord- 
wine  In  the  region  of  each  AIC  control  point. 
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The  integration  of  the  pressure  coefficient,  Eo.3.2  in  the  region 
of  the  AIC  point  leads  to  the  force 


q 


C  (I  ,  n  )  d£  dn 


(3.39) 


betting  '!  *  cos  4>  and  £  «  cos  0  In  Kqn.  3.2  arid  3.7  we  obtain 


Oa  4>0 
i 


f  e  q 


r*  «  ^  f 

?  )  )  n  /  j  si n^4>  s in  0  P  (  a )  f  (  £  )  d«j>  dO 

n  <  -  L-j  nm  I  |  m  '  n  '  * 


(  3.39A) 


where  0^,  and  0n  denote  the  chordwise  angular  measure  of  the  forward 


8f  4>j. 

denote  t 

and  aft  locations,  respectively,  of  the  pressure  region,  and  <i>.  and  <j>o 

denote  the  spanwise  angular  measure  of  the  inboard  and  outboard  end, 
respectively,  of  the  strip.  If  we  define  the  spanwise  integral 

Jfr  o 


sin  4>  Pm  (  jj.  )  d4> 


tJvV  ■ 

% i 

and  the  chordwise  integral 

-  _  ra 

Jn  (  ^f’  U  =  /  Bln0  fn  (  1  }  d( 

then  the  force  from  the  nm  mode  is 


(  3.3911) 


<  3.39C) 


urn 


q::  n  T  J 

nm  m  n 


(3.39D) 


and  in  matrix  form 

K  q  .*  /  |  1  J  j  \  a  j 

»  ,  |  m  n  |  |  nm  j 

so  that  the  Integration  matrix  j^Jis  £iyen  by 

i"l  -  i  (■x,.4„)'>n(  er,  l„) 


(3.39E) 


(3.39F) 
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The  integrals  are  easily  evaluated.  Consider  first  the  spanwise 
integration.  In  terms  of  the  spanwise  angular  coordinate  0  we  note 


P0  (  Jtt,  )  “1.0 


p„  <  «■ )  -  a  Vi  <  > 


*  cos  0  sin  m  0/sin  0  ,  1  <  m 


Then 


-4>c 


\  <*i>V  = 


sin  0  d0 


=  i  (cos"1Jlo  -  cos"1  ili  -Aq  1  -  jxo2  +  1-  R*  ) 


For  m  >  1 


I  = 
m 


_4*o 


sin  0  cob2  0  sin  m0  <30 


% 


from  which  4>n 

,■/ 


2  2 
sin  0  cos  0  d0 


3/2 


3/2 


-^(1-^!)  +  !„! 


nniJ  for  in  •!,  /  3, 


1 

m 


\  _ 

J  in  -.1 


_ L_ 

in  +1 


L 

m  -3 


.1 

m  +3 


aln  (m  -l)  0q  -  sin  (m  -l)  J 

J  sin  (m  +1)  0Q  -  sin  (m  +1)  0^ 

[sin  (m  -3)  0Q  -  sin  (m  -3)  <f>i  j 

l 

lain  (m  +3)  0Q  -  sin  (m  +3)  0..  j 


(3.40) 


(3.40A) 


(3.40B) 


(3. 40C) 


(  3.40C) 


((3.40D) 
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(3.40E) 


NOTE:  sin  n0  =  2  cos  0  sin  (n  -l)  0  -sin  (n  -2)  0 

so  that  sin  n0Q  =  2  sin  (n  -l)  0q  -  sin  (n  -2)  0q 

and  sin  n0^  =  2  sin  (n  -l)  0^  -  sin  (n  -2)  0^ 


Finally, 


f  9  5/2  9  5/2] 

3Ii-  |  bo (1  j 


3  L 

!  ,  , 

1  1  cos  n  -  1  cos  n .  -  1 

TT  I*°  u  _1  IT 


(3.40D) 
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/  1  -  TJ  2  -  21  Jl  - 

o  v  ~lo  Ji  y 


■^i 
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1 

5 


3/2 


DoU-Ho  )  -J*  (1-V ) 


3/2  1  ] 


Consider  next  the  chordwise  integral.  In  terms  of  the  chordwise 
angular  coordinate  0  we  note 


'o  <  *  >  55 


1  -  cos  0 
sin  0 


(3.41) 


fn  (  I  )  =  sin  n0,  n  21 


Jo  (  *  fl»  = 


(l  -  cos  0)  d  0 


i>r 


co r,  ^  -  cos 

n 


/ 


l-t*  ♦ 
■  a 


fTl ' 


(3.42) 


(3.43) 


For  n  '  .'] 


Jn  ’  2  (n  -1) 


fTPTTl) 


[  sin  (n  -1)  0a  -  sin  (n  -l)  0f  j 
|  sin  (n  +1)  0a  -  sin  (n  +l)  0f  j 


(3.44) 


find  for  n 


L 


till 

To  illustrate  the  format  of  I  B  J  consider  the  s  strip 


B_ 


^s^  Jn  *  ^2^ 


where  n  denotes  the  midpoint  of  the  s'1  strip,  and  1^,  I,  ...»  I  ^ 

denote  the  first  through  NC  chordvise  forces  on  the  strip.  Generalizing 
to  NS  strips  and  illustrating  the  dependence  on  n  and  a  we  have 


(3.46) 
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\  ( a-l>  Jo  < 

*0  <  *1>  Jo  <  «2> 
*o  <  V  Jo  <  V 

I0  (a2)d0  (fj.) 


*0  <  V  Jo  <  V 


Io  Jo  *  ^  NC^ 


(3.47) 


Ix  (iLg)  Jx  (  Ix) 
I1  ^2^  J1  ^  ^NC^ 


^jjS  JNS  ^  ^ 


Io  Jo  ^  ^  T1  (  nNS^  J1  ^1^ 


Z1  ^  J1  ^  ^NC^  XNS  ^  nNS^  JNS  ^  ? 


NC 
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DERIVATION  OF  SONIC  AND  SUPERSONIC  AERODYNAMIC  INFLUENCE  COEFFICIENTS 


The  supersonic  Mach  box  method  was  developed  by  Zartarian  and 

4 

Hsu  and  extended  to  intersecting  planar  lifting  surfaces  by  Moore  and 

5  6 

Andrew  .  The  sonic  box  method  was  developed  by  Rodemich  and  Andrew  , 

A  further  extension  to  the  wing-tail  combination  was  made  by  Moore 
and  Park  and  refined  by  Andrew3.  (NOTE:  These  references  are  out¬ 
lined  in  Part  I,  Section  A  for  the  sonic  case,  and  in  Part  VI,  Section  A 
for  the  supersonic  case. 

The  box  methods  lead  to  the  velocity  potentials  whose  strearrwise 
substantial  derivative  gives  the  pressure  coefficient.  The  solution  for 
the  velocity  potentials  may  be  written  £cf.  Eq.  (?)j 

W  ■  [aM2;1|  |w/vj  (3.48) 


The  pressure  coefficients  are  given  by 

|Cp|  =  (2/a  M)  |wj  |<i>|  (3.49) 

where  |wj  is  the  substantial  derivative  evaluated  at  the  box  centers. 

Since  the  boxes  are  all  small  the  aerodynamic  force  on  the  boxes  are 
approximately  given  by  the  product  of  the  box  pressure  and  its  area. 

The  forces  on  a  strip  (assumed  to  be  narrow)  are  regarded  as  acting  at 
the  strip  centerline  and  at  the  chordwise  centerline  of  each  spanwise 
line  of  boxes,  and  are  found  by  summing  spanwise  the  contribution  of 
each  box  to  the  strip. 


{Fa}  =  q  M  {Cp}  <3'50> 

The  elements  of  the  diagonal  area  matrix  j^j  consist  of  the  appropriate 
box  area  or  frnctiona  thereof  lying  on  the  strip  from  each  spanwise  line 
of  boxer.. 
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The  box  forces  are  converted  into  AIC  control  point  forces  by  the 
static  equivalence  of  the  box  forces  in  the  region  of  each  control  point, 
i.e.,  the  control  points  are  assumed  to  be  connected  by  a  series  of 
simple  beams,  hinged  at  each  control  point,  so  that  the  AIC  control  point 
forces  are  the  reactions  to  the  box  forces  distributed  along  the  series 
of  beams.  This  leads  to  approximate  generalized  forces  since  this 
method  of  representing  the  chordwise  deflections  is  not  consistent  with 
that  used  in  the  jV]  matrix.  BDwever,  it  leads  to  a  more  physically 
meaningful  distribution  of  AIC  forces,  and  the  resulting  approximation 
to  the  generalized  forces  in  an  arbitrary  vibration  mode  is  sufficiently 
accurate.  Denoting  the  statically  equivalent  transformation  matrix  by 
j^Tj  we  have  the  AIC  forces 


H  '  M  CFa} 

-  *  [*]  M  [°p] 


(3.51a) 

(3.51b) 


If  the  foregoing  equations  are  combined,  including  Eq.  3.8S  the 
AIC  control  forces  are  related  to  the  control  point  deflections  through 


jp)  ,  <*«/.. ">  [*]  [a]  [k„]  [a^J  [w]  (hj 

and  leads  to  the  AIC's  by  comparison  with  Eq.  3.1. 

[°h]  ■  W’ M  [T]  M  N  hi  >.i]  M 


(3.52) 


(3.53) 
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PART  IV 


SECTION  A 


TECHNICAL  DISCUSSION  OF  THE  SUBSONIC 
KERNEL  FUNCTION  METS5D 


INTRODUCTION 

In  defining  aerodynamic  influence  coefficients,  the  aerodynamic 
loads  are  first  derived  by  the  well  documented  "Kernel  function"  method, 
and  the  resulting  loads  are  then  converted  to  aerodynamic  influence 
coefficients.  The  procedure  followed  uses  results  from  any  of  the 
pertinent  papers,  and  extends  the  analysis  to  cover  the  problem  of 
tandem  surfaces. 

The  first  published  numerical  procedure  for  solving  the  subsonic 
pressure  distribution  problem  for  isolated  planar  lifting  surfaces  under¬ 
going  simple  harmonic  motion  was  developed  at  NASA's  Langley  Research 
Center  by  Watkins,  Runyan,  and  Woolston  (Reference  1).  Hsu  (Reference  7) 
significantly  advanced  the  logical  development  of  the  method  when  he 
established  an  optimum  set  of  collocation  and  integration  points. 

Rodemich  (Reference  8),  and  later  Landahl  (Reference  9),  have  presented 
expressions  for  the  kernel  that  are  very  much  simpler  than  those  previously 
used,  and  they  take  less  time  to  evaluate.  A  further  advance  was  made 
by  Rodden  and  Revell  (Reference  10),  who  described  the  matching  of  the 
boundary  conditions  with  the  least  squared  error.  Rowe  (Reference  ll) 
has  shown  that  to  obtain  sufficiently  accurate  results  using  Hsu's 
procedure,  an  extremely  large  number  of  collocation  points  must  be  used 
and  when  this  is  done,  using  the  spanwise  pressure  function  that  Hsu  and 
Watkins  used,  the  dowriwnsh  matrix  becomes  ill-conditioned.  Rowe  overcame 
this  problem  by  using  a  Fourier  series  for  the  spanwise  pressure  function. 

The  present  method  uses  Hsu's  net  of  collocation  points,  Rodemich's 
expression  for  the  kernel,  Rodden  nnd  Revell 's  idea  for  matching  the 
boundary  conditions,  and,  like  Rowe  does,  it  uses  an  orthogonal  set  of 
functions  for  the  npnnwloe  pressure  function.  Hbwever,  the  present  method 
utilizes  these  developments  in  a  way  that  has  not  been  previously  published. 
Tlic  techniques  employed  result  in  speeds  nnd  accuracy  not  previously 
al.l.a  i  ned . 
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TANDEM  SURFACE  VIBRATION  IN  SUBSONIC  FIOW 


The  method  presented  herein  vas  developed  for  application  to  a 
missile  with  a  very  low  aspect  ratio,  trapezoidal  wing  and  a  downstream 
rectangular  control  surface  lying  in  the  wake  of  the  wing. 

The  diameter  of  the  body  of  the  missile  is  considered  large  enough  to 
act  as  a  reflection  plane  for  acoustic  signals  that  emanate  from  points 
on  the  lifting  surfaces. 

The  problem  to  be  considered  is  that  of  determining  the  air  loads 
on  the  wing  and  control  surface  which  are  induced  by  a  simple  harmonic 
motion.  The  surrounding  fluid  is  assumed  to  be  compressible,  ir.viscid 
isentropic,  and  icrotational.  The  perturbation  potential  is  used,  and 
the  problem  is  further  linearized  by  applying  boundary  conditions  at 
the  mean  (s  *  0)  surface.  Thickness  effects  are  ignored.  With  these 
hypothesis,  it  is  well  known  that  an  integral  relation  exists  between 
the  pressure  discontinuity  over  the  surface  z  »  0  and  the  dovnwash 
over  the  same  plane.  If  the  downwash  and  pressure  difference  are 
representable  in  the  form 


Wx  (*,  y,  0,  t)  -  W  (x,  y)e  iwt 
AP1  (Sc,  y,  0,  t)  a  AP  (x,  y)e  iwt 


(4.1) 


The  integral  relation  becomes 


U 


TTu ~ 


^  (  §»  m,  x  - 1 ,  y-n)  d£  dn  (4.2) 


The  Integral  extends  over  the  plane  z  ■  0,  but  the  integrand  is 
zero  except  over  the  wing  and  control  surface.  The  kernel  function  K  is 
strongly  nlngular  and  integration  in  a  spanwise  direction  requires  the 
use  of  the  finite  part  concept.  This  is  indicated  by  the  cross  on  the 
Integration  sign.  Equation  (4.2)  then  represents  the  integral  equation 
of  the  nyatem  wherein,  given  W  (x,  y)  over  the  wing  and  control  surface, 
AP  muni  then  be  determined. 
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The  function  X  (k,  M,  x- §  ,  y-q  )  is  represented  here  in  the 

form 

K  (  M,  x-£  ,  y-n  )  *  e  U  (x  ~£  )  K-  i 

(y  -  n  )2  (4.3) 

where 


M  =  Mach  number 

p2  =  1  -  M2 

K1  ^1^  is  modified  Bessel  function  of  the  second  kind  of  order 
one  and  argument,  k^. 

Thin  I'orm  for  the  kernel  function  has  been  used  in  References  8  and 
9.  While  It  la  not  identical  to  that  of  Reference  1,  it  may  be  obtained 
directly  from  that  equation  by  substitution  of  the  integral  representation 
of  the  modified  Bessel  and  Struve  functions.  The  second  order  singularity 
in  the  npanwlso  variable  in  Equation (4.3)  requires  the  use  of  the  "finite 
part"  technique  of  Hadamard. 
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PRESSURE  DISTRIBUTION 

As  in  many  earlier  studies,  the  pressure  distribution  is  approxi¬ 
mated  as  the  sum  of  a  series  of  functions  which  have  the  proper  behavior 
as  inferred  from  steady  state  and  two-dimensional  solutions.  This  behavior 
includes  a  Kutta  condition  at  the  trailing  edge,  a  square  root  singularity 
at  the  leading  edge,  and  a  half  ordered  zero  at  the  top  for  each  surface. 
The  pressure  distribution  on  each  surface  is  then  approximated  in  the  form 


ap(g>  n) 
\  pu2 


* 


v's2  -  Tp 


N  M 

fn  ® 


n=0  n=Q 


(4.5) 


where 
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The  functions  are  Chebyshev  Polynomials  and  are  introduced  for 
purposes  of  convenience. 

For  a  trapezoidal  wing,  the  lines  f  »  constant  and  JJ»  constant  are 
not  orthogonal.  This  transform tion  maps  each  surface  into  a  square  in  its 
(  1.  H  )  plane. 


The  fundamental  integral  Equation  (4.2)  is  inverted  by  substituting 
Equation  (4.5)  into  Equation  (4.1),  and  determining  the  coefficient  a 

run 

so  that  Equation  (4.2)  is  equilibrated  at  a  designated  set  of  collocation 
points.  This  equilibrium  is  represented  schematically  in  the  form 


(4.2) 


The  left  hand  side  of  Equation  (47)  represents  the  prescribed 

downwash  at  the  chosen  set  of  collocation  points.  D  a  is  the  effect 

run  nm 

on  downwash  of  the  corresponding  term  in  the  pressure  series  expansion. 

The  superscripts  W  and  C  designate  wing  and  control  surface,  respectively. 
The  left  superscript  on  the  D's  indicate  the  surface  on  which  the 
col  local. ion  point  is  located,  the  right  superscript,  the  surface  over 
which  the  Integral  ia  taken.  Then  (x,  y)  is  given  by 

2  f  1  - ?  f 1  K  (x-g,  y-n) 

•  '  fc  f  V1  -a  P«(3>  fn<f>  ~ — X”  ^‘8) 

J. i  J.i  (y  -n) 
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The  integrals  involved  in  Equation (4. g)  ate  evaluated  by  following 
the  sethods  of  Hsu,  So  this  procedure,  the  Cauas4iehler  quadrature  is 
used,  and  the  collocation  point*  are  selected  by  a  asihod  which  is  analogous 
to  using  the  Gauss-Hebter  quadrature  on  the  inverse  problem. 

In  brief,  this  technique  is  concerned  with  a  numerical  evaluation 
of  an  Integral  in  the  form 

i” 

Vfx)  £(x)  ix  «  ^  H,  £(x A 

J  r«i  J  i 

Where  W  (*.)  is  a  given  weighting  function,  f(x)  an  arbitrary 
function,  and  Hj  are  weighting  number*.  The  first  choice  of  points  X^. 
is  taken  to  be  chat  corresponding  to  which  the  approximate  integration 
would  be  exact  for  a  polynomial  of  degree  leas  than  or  equal  to  2n  -1. 

The  Integration  points  and  weighting  numbers  are  listed  in  Rope!  (Reference 
12)  page  2S3.  They  are;  for 
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where  .1  »  1,  n 
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The  inner  integral  of  Equation  (4.%  is  then  evaluated  by  using 
the  third  of  Equation  (4.9),  With  the  previously  defined  notation  then 


r«®*  ■  2  vA> 


(4.10) 


where 


s  _  /2k  •  1  v 

h  *  *coa  (gm  ") 


Hsu  pointed  out  that  if  the  collocation  points  are  interdigitated  accord¬ 
ing  to  the  formula 


*i  *  *  <srfi n) 


the  inner  integral  is  evaluated  with  the  least  squared  error  by  Equation 
(4.10). 


This  minimization  was  of  course  achieved  for  a  single  surface,  but 
the  convention  was  retained  for  the  tandem  surface  model. 
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Hsu  did  not  take  advantage  of  the  fact  that  the  number  of 
collocation  points  does  not  need  to  be  as  great  as  the  number  of  integratiori 
points.  In  the  present  method  the  computer  program  user  may  specify  the 
number  of  integration  points  to  be  any  positive  integer  times  the  number 
of  chordwise  collocation  points,  so  long  as  the  number  of  integration 
points  is  40  or  less.  Sixteen,  or  more,  chordwise  integration  points 
are  recommended.  Accurate  results  can  be  assured  only  if  the  user 
specified  a  sufficient  number  of  quadrature  points  for  the  integral. 

Substitution  of  Equation  (4.10)  into  Equation  (4.8)  gives 


r,  “S 

Dnm"BT 


(x,  y  -  Tl)  dfl_ 


(4.11, 


where 


nn 


2L  +  T  pm^)  £  0-  -  sk2)  ua„i^  V*  -  Sk>  y  -  11);  l  s  n 

k«l 


The  integral  for  is  singular.  The  Gauss-Mehler  technique  is 

again  followed.  This  time  the  first  of  Equations  (4. 9 )  is  used  to  perform 

the  integration,  and  the  finite  part  modification  is  introduced.  With 

these  constraints,  Hsu  has  shown  that  the  evaluation  for  D  becomes 
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I 


and 


X 

°>  -  -a /  "  5)a? 

When  the  collocation  point  lies  on  the  control  surface  this  integral 
is  taken  over  the  wing,  from  the  leading  edge  (  f  =  -1)  to  the  trailing 
edge  C  |  =  +1),  and  over  the  control  surface,  from  the  leading  edge  to  the 
collocation  point  (  £;  -  x) .  The  integrand  of  this  expression  does  not 
vary  widely  over  the  region  of  integration  and  is  evaluated  with  sufficient 
accuracy  by  a  six-point  quadrature  formula.  It  has  been  made  a  part  of 
the  computer  program  and  may  not  be  controlled  by  the  user. 
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PART  IV  -  SECTION  B 

SUBSONIC  AIC  COMPUTER  PROGRAM  DESCRIPTION 


A  FORTRAN  IV  computer  program  is  presented  which  computes  subsonic 
unsteady  aerodynamic  influence  coefficients  for  a  variety  of  tandem  coplanar 
lifting  surface  configurations.  The  computer  solution  is  based  on  a  kernel 
function  formulation  which  satisfies  the  linearized  equations  of  motion  of  an 
inviscid,  isentropic,  compressible  and  irrotational  fluid.  The  analysis  is 
extended  to  include  interaction  effects  between  tandem  surfaces  and  wake  effects 
on  the  trailing  surface. 

The  various  configurations  which  can  be  analyzed  are  shown  in  Figure  4.1. 

The  vehicle  body  is  considered  as  a  reflection  plane  for  acoustic  signals 
emanating  from  the  aerodynamic  surfaces  thereby  giving  the  surfaces  a  plane 
of  symmetry,  parallel  to  the  free-stream  flow.  The  upstream  surface  (wing) 
must  have  an  unswept  trailing  edge  and  the  rectangular  trailing  surface  must 
have  the  same  spanwise  dimension  as  the  trailing  edge  of  the  wing.  The  down¬ 
stream  surface  lies  in  the  wake  of  the  wing;  any  non-negative  value  may  be  used 
for  the  gap  dimension.  A  single  surface; cannot  be  analyzed,  however,  an  option 
is  provided  which  eliminates  interaction  effects.  Thus  it  is  possible  to  generate 
AIC's  for  individual  surfaces  isolated  from  disturbances  in  the  flow  field. 

The  program  allows  up  to  40  AIC  control  points,  20  per  surface.  The  AIC 
stations  must  satisfy  the  following  requirements: 

(1)  Both  surfaces  must  have  the  same  number  of  spanwise  rows  of  control 
points.  The  chordwise  location  of  respective  rows  on  the  surfaces  need 
not  be  the  same. 

(2)  The  chordwise  rows  must  be  parallel  to  the  flow  stream. 

(3)  The  chordwise  rows  on  a  surface  must  have  the  same  number  of  control  points. 

(4)  The  control  points  in  each  spanwise  row  must  have  the  same  fractional 
chord  location. 

(5)  The  minimum  number  of  chordwise  or  spanwise  control  points  for  a  surface 
Is  two  and  the  maximum  number  is  ten. 

(6)  The  origin  for  the  AIC  station  coordinates  and  the  geometric  coordinates 
for  the  plunform  must  be  at  the  leading  edge  root  of  the  wing. 


I 

i 


Examples  of  acceptable  AIC  control  point  patterns  for  the  subsonic  program 
i  1 1  list  raLeil  in  Figure  4.2. 


are 


fRAPSZOIOAl 


VVW  I 


FIGURE  k.l 


TANDEM  COPLANAR  CONFIGURATIONS  AT  SUBSONIC  MACH  NUMBER 


The  computer  program  consists  of  a  main  program  (MAIN)  and  24  subroutines 
and  functions  subprograms.  Execution  begins  with  MAIN  calling  KFDA  which  reads 
input  data  and  stores  this  information  in  core.  The  program  then  calls  TRAMP 

.  The  £wj  matrix  serves 
two  functions.  It  relates  the  collocation  stations  of  the  unsteady  aerodynamics 
to  the  control  stations  of  the  AIC  matrix  and  it  serves  as  a  substantial 
derivative  operator.  Subroutines  called  by  TRAMP  are  CMAT,  SMAT,  TMAT,  BMAT, 

RMAT  and  MINV. 

The  next  item  computed  is  the  kernel  function  matrix,  j^j.  The  subroutine 
CORD  is  called  for  each  unsteady  aerodynamic  collocation  station  and  constructs 
the  kernel  function  matrix  which  is  dependent  only  on  the  relative  location  of 
the  collocation  stations  and  the  Mach  number-reduced  frequency  combination. 

The  pressure  coefficients  are  found  from  the  relation 

M =  M'1  M 

The  program  has  been  written  such  that  the  number  of  spanwise  and  chordwise 
pressure  coefficient  terms  and  the  number  of  spanwise  and  chordwise  collocation 
stations  for  the  unsteady  aerodynamics  matches  the  respective  number  of  AIC 
control  points  on  each  surface.  Thus  the  Kernel  function  matrix  is  square  and 
its  inverse  is  computed  directly.  The  subroutines  employed  for  this  operation 
are  CGRED  and  XLSQ. 

After  deriving  the  pressure  coefficients,  the  pressure  terms  are  integrated 
spanwise  and  chordwise  to  obtain  the  force  acting  at  each  AIC  control  station. 

The  resulting  matrix,  after  it  is  multiplied  by  a  non-dimensionalizing  factor, 
is  the  final  AIC  matrix,  'jc^J .  By  de££alti'«int!te’ AIC>  matrix  relates  forces  to 
displacements  through  the  equation 

{f}  -  Pi.2  •  [chJ 

The  semi- chord  of  the  wing  root  is  used  hs  the  reference  chord,  br.  s  is  the 
semi-span  of  the  wing  (and  tail)  and  co  is  the  oscillatory  frequency  in  radians/ 
sec.  This  final  phase  of  the  subsonic  AIC  development  uses  the  subroutines  AICS, 
KORGK,  AKCCOS,  MINTS,  and  MINTC. 


which  generates  the  substantial  derivative  matrix, 
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1.0 


PROCESSING  REQUIREMENTS 


The  input  and  output  files  used  by  the  program  are  05  and  06,  respectively. 
All  read  and  write  statements  are  contained  in  the  main  program  (MAIN)  and  the 
subroutines  KFDA  and  KOUT.  Peripheral  tape  and  disc  units  are  not  used  by  the 
program.  Approximately  40,000  cells  of  core  storage  is  required.  A  standard 
input  form  of  six  12-column  fields  per  card  is  used  by  the  program.  Floating 
point  numbers  (6E12.5  format)  may  lie  anywhere  within  the  appropriate  field, 
but  fixed  point  numbers  (6112  format)  must  be  right  adjusted.  Detailed 
instructions  for  data  input  are  given  and  listings  of  data  for  sample  problems 
are  provided. 


2.0  INPUT  INSTRUCTIONS 


Instructions  for  preparing  input  data  for  the  subsonic  AIC  computer  program 
are  presented  here.  The  field  location  and  format  for  each  input  quantity  is 
specified.  Any  set  of  units  may  be  used  for  geometric  dimensions  and  acoustic 
velocity  as  long  as  they  are  consistent  e.g.,  if  feet  is  used  for  length  then 
the  acoustic  velocity  must  have  dimensions  of  feet  per  second. 

1.  Streamwise  Coordinates  (6E12.5  format) 


Column 

Name 

Item 


(1)  X(l)  Wing  root  leading  edge  streamwise  coordinate 

(2)  X(2)  Wing  tip  leading  edge  streamwise  coordinate 

(3)  X(3)  Wing  trailing  edge  streamwise  coordinate 

(4)  X(4)  Control  surface  leading  edge  streamwise  coordinate 

(5)  X(5)  Control  surface  trailing  edge  streamwise  coordinate 

The  technique  for  generating  various  configurations  is  shown  in  Table  4.1. 

The  origin  for  the  planfcrm  and  AIC  station  coordinates  must  be  at  the 
leading  edge  root  of  the  wing  therefore  X(l)  and  Y(i),  described  below, 
must  always  be  0.0. 

2.  Spanwise  Coordinates  and  Acoustic  Velocity  (6E12.5  format) 

Column 
Name 
Item 


(1)  Y(l)  Wing  root  spanwise  coordinate 

(2)  Y(2)  Wing  leading  edge  spanwise  coordinate 

(3)  Y(3)  Wing  (and  control  surface)  tip  spanwise  coordinate 

(4)  SOUND  Acoustic  velocity  for  altitude  at  which  analysis  is  performed 


1  1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

CBM 

mm 

SOUND 

IBM 

fSMI 

wsm 

CBM 

1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

BBC! 

MM 

Hi 

mm 

wsm 

IBM 

wsm 

CBM 

CBM 

TABLE  4.1  -  OPTIONAL  CONFIGURATIONS 


l 


Configuration 

Chordwise  Coordinates 

Spanwise  Coordinates 

Rectangular 

X(l)  =  0.0 

Y(l)  =  0.0 

X(2)  =  0.0 

Y(2)  -  0.0 

X(3)  >  0.0 

X(4)  >  X(3) 

X ( 5)  >  X(4) 

Y(3)  >  0.0 

Delta 

X(l)  =  0.0 

Y(l)  ■»  0.0 

X(2)  >  0.0 

Y(2)  =  0.0 

X(3)  =  X(2) 

X(4)  >  X(3) 

X(5)  >  X(4) 

Y(3)  >  0.0 

Trapezoidal 

X(l)  =  0.0 

Y(l)  =  0.0 

X(2)  >  0.0 

Y(2)  >  0.0 

X(3)  =  X(2) 

X(4)  >  X(3) 

X(5)  >  X(4) 

Y(3)  >  Y(2) 

Trapezoidal 

X(l)  =  0.0 

Y(l)  =  0.0 

(Cropped) 

X(2)  >  X(l) 

Y(2)  >  0.0 

X(3)  >  X(2) 

X(4)  >  X(3) 

X(5)  >  X(4) 

Y  (3)  >  Y(2) 

Delta 

X(l)  -  0.0 

Y(l)  =  0.0 

(Cropped) 

X(2)  >  0.0 

Y(2)  =  0.0 

X(3)  >  X(2) 

X(4)  >  X(3) 

X(5)  >  X(4) 

Y (3)  >Y(2) 
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3.  General  Information  (6112  format) 
Column 
Name 
Item 


(1)  NMACH  Number  of  Mach  numbers  (max.  6) 

(2)  KF  Option  to  input  either  frequencies  or  reduced  frequencies: 

KF  ®  0  frequencies 

KF  =  1  reduced  frequencies 

(3)  NFREQ  Number  of  frequencies  or  reduced  frequencies  at  each  Mach 

number  (max.  10) 

(4)  LCOLL  Option  to  print  aerodynamic  and  AIC  collocation  station 

coordinates: 

LCOLL  =  0  do  not  print 

LCOLL  =  1  print  collocation  station  coordinates 

(5)  LPUNCH  Option  to  punch  AIC  matrix  on  cards: 

LPUNCH  =  0  do  not  punch 
LPUNCH  *  1  punch  AIC  matrix  for  wing  only 
LPUNCH  =  2  punch  A.IC  matrix  for  control  surface  only 
LPUNCH  *  3  punch  individual  AIC  matrix  for  wing  and 
control  surface 

LPUNCH  =  4  punch  total  AIC  matrix  of  combined  wing  and 
control  surface 

The  AIC  matrix  is  punched  by  rows  with  a  1P6E12.5  format.  Each  row  of 
the  matrix  begins  on  a  new  card. 


4.  General  Information  (6112  format) 


Column 

Name 

Item 


1-12 

13-24 

25-36 

17-48 

49-60 

61-72 

NWCX 

NCCX 

NIONCX 

NIY 

ISOLAT 

■m 

mm 

ISHMI 

IBM 

0) 

(2) 

(3) 


(4) 

(5) 


NWCX  Number  of  chordwiae  collocation  stations  on  wing  (2  5NWCXS  10) 

NCCX  Number  of  chordwiae  collocation  stations  on  control  surface 

(2  £  NCCX  £  10) 

NIONCX  Factor  for  number  of  chordwise  integration  stations.  Set 
NIONCX  such  that  NI0NCX*NWCX  and  NI0NCX*NCCX  are  greater 
than  15  (but  less  than  40)  to  insure  sufficient  quadrature 
points  for  accurate  integration  of  the  Kernel  function.  If 
AIC's  are  desired  for,  say,  the  wing  only,  set  NIONCX  such 
that  NI0NCX*NWCX  is  greater  than  15  and  set  NCCX  equal  to  2 
to  minimize  computing  time. 

NIY  Number  of  spanwise  collocation  stations  (wing  and  control 

surface) 

ISOLAT  Option  to  isolate  wing  and  control  surface: 

ISOLAT  **  0  interference  and  gap  effects  considered 
ISOLAT  ■  1  surfaces  are  isolated.  AIC's  will  be  for 

individual  surfaces  with  no  coupling  effects. 


5.  Mach  Numbers  (6E12.5  format) 


Column 

Name 

Item 


1-12 

13-24 

25-36 

37-48 

49-60 

61-72  | 

lliilBI 

SMB 

HBRS 

mm 

KH 

mm 

Hi 

hi 

I5H 

(1)  FMACH  (1)  Mach  number 

(2)  FMACH  (2)  Mach  number 


(NMACH)  FMACH (NMACH)  Mach  number 

NMACU  values  of  Mach  number  must  be  input  (see  Part  3,  Item  1).  Mach 
numbers  must  be  greater  than  zero  and  less  than  0.95. 
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6.  Frequencies  or  Reduced  Frequencies  (6E12.5  format) 


I 

1 


Column 
Name 
Item 

If  KF=0, input  NFREQ  values  of  frequency  (cps) .  If  KF=1,  input 
NFREQ  values  of  reduced  frequency  (k^  =  twby/U  where  =  semi-chord  of  wing 
root,  U  »*  free  stream  velocity,  and  w .» oscillatory  angular  frequency  in 
radians/sec).  Frequencies  (and  reduced  frequencies)  may  not  be  zero. 

(1)  1REQ(1)  f(cps)  or  k_ 

i. 

(2)  FREQ (2)  " 


"1-12 

13-24 

25-36 

'37-48 

49-60 

61-72 

FREQ(l) 

FREQ (2) 

FREQ (3) 

FREQ (4) 

FREQ (5) 

FREQ (6) 

jn _ 

in _ 

sn _ 

_ 

11} _ 

i& _ 

(NFREQ)  FREQ (NFREQ) 

For  NFREQ  >6,  continue  input  on  new  card. 


7.  Spanwise  Location  of  AIC  Stations  on  Wing  (6E12.5  format) 

Column 
Name 
Item 

(1)  YAIC(1,W)  Spanwise  coordinate  of  first  row  of  AIC  collocation 

stations  on  wing. 

(2)  YAIC(2,W)  Spanwise  coordinate  of  second  row  of  AIC  collocation 

stations  on  wing. 


1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

YAIC(l.W) 

YAIC(2,W) 

YAIC(3,W) 

YAIC (4, W) 

YAIC(5,W) 

YAIC  (6  ,W) 

jn _ 

ilL. 

in _ 

_ 

in _ 

(6) 

(NIY)  YAIC(NIY,W)  Spanwise  coordinate  of  last  row  of  AIC  collocation 

stations  on  wing 

CoLlocation  station  rows  are  numbered  from  root  to  tip.  For  NIY  >  6, 
continue  Input  of  YAIC  (7,W)  to  YAIC(NIY,W)  on  new  card(s). 
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Spanwise  Location  of  AIC  Stations  on  Control  Surface  (6E12.5  format) 


Column 

Name 

Item 


(1)  YAIC(1,CS)  Spanwise  coordinate  of  first  row  of  AIC  collocation 

stations  on  control  surface. 

(2)  YAIC(2,CS)  . . .  second . 


1  1-12 

13-24 

25-36 

EBB H 

49-60 

61-72 i 

EfiHBS 

tmm 

ISflHH 

ffBSHHK 

INH 

(SHH 

rrmm 

(NIY)  YAIC(NIY,CS)  .  last . 

Collocation  station  rows  are  numbered  from  root  to  tip.  If  NIY>6, 
continue  input  of  YAIC(7,CS)  to  YAIC(NIY,CS)  on  new  card(s). 

9.  Streamwise  Location  of  AIC  Stations  on  Wing  (6E12.5  format) 

Column 
Name 
Item 


(1)  XAIC(W,1,1)  Streamwise  coordinate  of  first  AIC  collocation  station 

in  first  row  on  wing. 

(2)  XAIC(W,  1,2)  . . .  second . 


1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

1S0BSHBB 

EfflHHEl 

B19I9BBI 

•  •  • 

»  •  • 

•  •  • 

jn _ _ 

12) _ 

jzl _ 1 

181 

ft  1 

rm 

(NXWING*NYWING) 

XAIC(W,NYWING,NXWING) .  last . last  row . 

Si reanvwisc*  numbering  sequence  if  from  leading  edge  to  trailing  (see 
Figure  4.2).  Continue  input  of  values  for  each  row  immediately  after 
the  last  value  of  the  preceding  row,  do  not  begin  input  of  each  new 
row  on  new  card. 


10.  Chordwlse  Location  of  AIC  Stations  on  Control  Surface  (6E12.5  format) 

Column 
Name 
I  lorn 

Procedure  to  Input  streamwise  coordinate  location  of  AIC  stations 
on  control  surface  Is  the  same  as  wing. 


13-24 

25-36 

37-48 

^49-60 

61-72 

IBRSHffl 

•  *  • 

. . . 

•  •  • 

(!) 

ftHB 

Ua _ 

1 

IqI 
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3.0  Sample  Problems 

Three  sample  problems  are  given  to  demonstrate  the  operation  of  the 
subsonic  AIC  program  and  to  illustrate  the  options  available.  Trapezoidal, 
cropped  trapezoidal,  and  delta  configurations  are  analyzed.  Explanation  of 
input  parameters  and  computer  listings  of  input  data  cards  and  computer  output 
arc  given  for  each  problem. 

Sample  Problem  1 

A  cropped  Delta  wing-rectangular  control  surface  combination  is  analyzed  for 
M  =  0.5  and  k^  =  0,10,  The  planfortn  geometry  and  AIC  control  station  locations 
are  shown  in  Figure  4.4.  The  dimensional  unit  used  for  length  is  feet  for  this, 
particular  case,  therefore  the  acoustic  velocity  is  entered  with  units  of  feet 
per  second.  NIONCX  is  set  equal  to  8  which  makes  NIONCX*NWCX  and  NIONCX*NCCX 
greater  than  15,  thereby  insuring  sufficient  chordwise  quadrature  points  for 
accurate  chordwise  integration  of  the  kernel  function. 

The  surfaces  have  5  spanwise  collocation  stations  and  the  wing  has  3 
chordwise  stations  while  the  control  surface  has  2  chordwise  stations. 
Summarized  below  are  input  parameters.  A  listing  of  the  data  input  cards  and 
computer  output  follows. 


X(l)  =  0.0' 

Y(1 )  =  0.0' 
SOUND  =  1116.87 
NMACH  =  1 
KF  -  1 
NFREQ  =  1 
LCOLL  =  1 
LPUNCH  =  4 
NWCX  «  3 
NCCX  ==  2 

NIONCX  =  8 

NIY  •*  5 
I  SOT  AT  -  0 
FMACII(I)  ■  0.5 
KKKQ(I)  -  0.1 


X(2)  =  1.0'  X(3)  =  2.0'  X(4)  »  3.0'  X(5)  =  4.0' 

Y(2)  =  0.0'  Y(3)  =  2.0' 

ft/sec  (analysis  for  Sea  level) 

number  of  Mach  numbers 

Input  reduced  frequency 
Number  of  reduced  frequencies 
Print  collocation  station  coordinates 
Punch  total  AIC  matrix  on  cards 

Number  of  chordwise  AIC  collocation  stations  on  wing 

Number  of  chordwise  AIC  collocation  stations  on 
control  surface 

Factor  for  determining  number  of  chordwise  integration 
stations 

Number  of  spanwise  AIC  collocation  stations 
Surfaces  are  not  isolated 
Mach  number 
Reduced  frequency 


-48- 


YAIC(1,W)  =  0 
YAIC(4,W)  *  1. 

YAIC(1,CS)  « 
YAIC(4,CS)  = 

XAIC ( 1 , 1 ,  W)  - 
XAIC (2, 1,W)  = 
XAIC(3,1,W)  = 
XAIC(4,1,W)  = 
XAIC(5,1,W)  = 

XAIC(1 , 1 ,CS)  = 
XAIC(2, 1 ,CS)  = 
XAIC(3,1,CS)  = 
XAIC(4,1,CS)  - 
XAIC(5, 1 ,CS)  = 


2' 

YAIC(2,W)  =  0.6' 

YAIC(3,W)  r  1.0’ 

4' 

YAIC(5,W)  =1.8' 

1.1' 

YAIC(2,CS)  =  0.5' 

YAIC(3,CS)  =  1.0' 

..5' 

YAIC(5,CS)  =  1.9' 

0.575’ 

XAIC ( 1 , 2 , W)  -  1.050' 

XAIC(1 ,3,W)  =  1.525' 

0.725' 

XAIC(2,2,W)  =  1.150' 

XAIC(2,3,W)  =  1.575' 

0,875' 

XAIC(3, 2,W)  =  1.250' 

XAIC(3,3,W)  =  1.625' 

1.025' 

XAIC (4, 2, W)  =  1.350' 

XAIC(4,3,W)  =  1.675' 

1.175' 

XAIC ( 5 , 2 ,  W)  =  1.450' 

XAIC(5,3,W)  =  1.725' 

3.25' 

XAIC(1,2,CS)  =  3.75' 

3.25' 

XAIC(2,2,CS)  =  3.75' 

3.25' 

XAIC(3,2,CS)  =  3.75' 

3.25' 

XAIC(4,2,CS)  «  3.75' 

3.25' 

XAIC (5, 2, CS)  =  3.75' 
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f)  A  TA  HARD  COLUMN  NUMBER 
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Sample  Problem  2 

A  cropped  trapezoidal  wing  is  analyzed  for  M  ®  0.5  and  *  0.10.  Since 
the  subsonic  AIC  program  requires  a  control  surface,  a  dummy  surface  is  added 
with  a  minimal  number  of  chordwise  collocation  stations  (2)  to  minimize 
computing  time.  The  surfaces  are  isolated  and  the  option  LPUNCH  is  input  as 
1  to  punch  the  AIC  matrix  for  the  wing  only.  The  wing  has  4  chordwise  and  4 
spanwise  collocation  stations.  Planform  geometry  and  AIC  control  station 
locations  are  shown  in  Figure  J+.6.  Summarized  below  are  input  parameters.  A 
listing  of  the  data  input  cards  and  computer  output  follows. 


X(l)  =  0.0’  X(2)  =  1.0*  X(3)  =  2.0*  X(4)  *  3.0'  X(5)  *  4.0' 

Y(l)  =  0.0'  Y(2)  =  1.0'  Y(3)  =  2.0' 


SOUND  =  1116.87  ft/sec 

NMACH  =  1 

KF  *  1 

NFREQ  =  1 

LCOLL  =  1 

LPUNCH  =  1 

NWCX  =  4 

NCCX  =  2 

NIONCX  =  4 

NIY  =  4 
ISOLAT  =  1 
FMACH  (1)  =0.5 
FREQ(l)  =  0,10 


Number  of  Mach  numbers 
Input  reduced  frequency 
Number  of  reduced  frequencies 
Print  collocation  station  coordinates 
Punch  AIC  matrix  for  wing 

Number  of  chordwise  AIC  collocation  stations  on  wing 

Number  of  chordwise  AIC  collocation  stations  on 
control  surface 

Factor  for  determining  number  of  chordwise  integration 
stations 

Number  of  spanwise  AIC  collocation  stations 
Isolate  wing  and  control  surface 
Mach  number 
Reduced  frequency 


YAIC(1 ,W)  -  0.2'  YAIC(2,W)  -  0.7'  YAIC(3,W)  =  1,3' 

YAIC(4,W)  -1.8' 


YAIC(t.CS)  -  0.2’  YAIC(2,CS)  *  0.7'  YAIC(3,CS)  =  1.3' 

YAtC(4,CS)  ~  1.8' 
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xaic(i,i,w)  =  o.io' 

XAIC(1,2,W)  »  0.70' 

XAIC(1,3,W)  =  1.30’ 

XAIC(1,4,W)  *  1.90* 

mc(2,i,«)  -  o.io' 

XAIC(2,2,W)  »  0.70' 

XAIC(2,3,W)  -  1.30' 

XAIC(2,4,W)  «  1.90' 

XAIC(3,1,W)  =  0.38* 

XAIC(3,2,W)  «*  0.90’ 

XAIC(3,3,W)  =  1.405' 

XAIC(3,4,W)  =*  lo915' 

XAIC(4,1,W)  «  0.86' 

XAIC(4,2,W)  -  1.22' 

XAIC(4,3,W)  =  1.56' 

XAIC(4,4,W)  =*  1.94' 

XAIC(1,1,CS)  =  3.25' 

XAIC(1,2,CS)  «  3.75' 

XAIC(2,1»CS)  =  3.25’ 

XAIC(2,2,CS)  »  3.75* 

XAIC(3,1,CS)  =  3.25' 

XAXC(3,2»CS)  -  3.75' 

XAIC(4,1,CS)  =  3.25' 

XAIC(4,2,CS)  *  3.75' 
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SUBSOMC  SAMPLE  PROBLEM  *Z 


DATA  CARD  COLUMN  NUMBFR 
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A  45°  delta  wing-control  surface  combination  is  analyzed  for  M  *  0.5  and 
f  »  5.0  cps.  There  are  4  spanwise  and  4  chordwise  AIC  control  stations  for 
both  the  wing  and  control  surface.  The  planform  geometry  and  AIC  control 
station  locations  are  shown  in  Figure ^.8.  The  input  parameters  are  suianarized 
below  and  a  listing  of  the  data  input  cards  and  computer  output  follows. 


X(l)  =  0.0'  X(2)  *  2. 0‘  X(3)  *  2.0'  X(4)  =  3.0'  X(5)  =  4.0» 

Y(l)  =  0.0'  Y(2)  *  0.0'  Y(3)  =  2.0' 


SOUND  «  1116.87  ft/sec 

NMACH  =  1 

KF  =  0 

NFREQ  »  1 

LC0LL  =  1 

LPUNCH  =  4 

NWCX  *  4 

NCCX  =  4 

NIONCX  =  4 

NIY  =  4 
ISOLAT  =  0 
FMACH(l)  =  0.5 
FREQ(l)  =  5.0 


Number  of  Mach  numbers 

Input  frequencies 

Number  of  frequencies 

Print  collocation  station  coordinates 

Punch  cards  for  total  AIC  matrix 

Number  of  chordwise  AIC  collocation  stations  on  wing 

Number  of  chordwise  AIC  collocation  stations  on 
control  surface 

Factor  for  determining  number  of  chordwise  integration 
stations 

Number  of  spanwise  AIC  collocation  stations 
Surfaces  are  not  isolated  , 

Mach  number 
Frequency 


YAICU.W)  =  0.2'  YAIC(2,W)  =«  0.6'  YAIC(3,W)  -  1.0' 

YAIC(4,W)  «  1.4' 


YAIC(1,CS)  -  0.25' 
YAIC(4,CS)  «  1.75' 


YAIC(2,CS)  -  0,75'  YAIC(3,CS)  -  1.25' 


XAIC(1, 1,W)  -  0.56* 
XAIC(1,4,W)  «  1.64’ 
XAIC(2,1,W)  »  0.88* 
XAIC(2,4,W)  »  1.72' 
XAIC(3,1,W)  *  1.20' 
XAIC(3,4,W)  =1.80' 
XAIC<4,1,W)  =  1.52' 
XAIC(4,4,W)  «  1.88' 

XAIC(1,1,CS)  =  3.125* 
XAIC(1,4,CS)  «  3.875* 
XAIC(2, 1,CS)  =  3.125* 
XAIC(2,4,CS)  =  3,875* 
XAIC(3,1,CS)  =  3.125* 
XAIC(3,4,CS)  =  3.875’ 
XAIC(4,1,CS)  »  3.125* 
XAIC(4,4,CS)  =  3.875* 


XAIC(1,2,W)  =  0.92* 
XAIC(2t2,W)  *  1.16* 
XAIC(3,2,W)  -  1.40* 
XAXC(4,2,K)  =  1.64* 

XAIC{1,2,CS)  -  3.375’ 
XAXC(2,2,CS)  -  3.375’ 
XAIC(3,2,CS)  -  3.375’ 
XAIC(4,2,CS)  -  3.375’ 


XAIC(1,3,W)  =  1.28’ 
XAIC(2,3,W)  =  1.44’ 
XAIC(3,3,W)  =  1.60* 
XAIC(4,3,W)  =  1.76' 

XAIC(1,3,CS)  -  3.625’ 
XAIC(2,3,CS)  =  3.625’ 
XAIC(3.3,CS)  =  3.625* 
XAIC(4,3,CS)  =  3.625* 


DATA  CARD  COLUMN  NUHUFR 
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PART  IV  -  SECTION  B4.0 


LISTING  OF  SUBSONIC  AIC  COMPUTER  PROGRAM 


CMA  IN 


HA  JN 

COMPLEX  A* A A . ANH.CZtRO.WASH.  A  IC 

51  f  Hf  NS  I  ON  Tf  MPMO  ,4ii  ),  AfCMIl  .O'1 ) .  HASHMli )  ,FM(4n,40  ) 

COHHON/Cl  /  A<  «!»  ),AA(An,no),  ANH(  4  o  .  4H  )  .CZfcRO 

COMMON/C?/HKFR(  ->0  )  .  ZKFR(?U>*FMAC'H(^)»FRfcO(10  )  , NOH(b  )  .  IL  (‘>o  )* 
l  HCOR(ft)  .ZrOR(6).HXnMN(i  1  )  *  UBCNf  J 1  ).UBIN(in.H!(90). 

?  XM*»>»  Ye(.O.UX(lO).UY(iii)»WXlHNni  ), SIXMO  /* )  #  SCX  ( 1  U  •  2  )  * 

.5  E  T  A  (  1  ) 

COMMON/C-V  Y<  i  l».XAin(ln#t».i).YAFC<lfl#?).B(4fl#4«),R(‘Wi.Ail). 

J  C(»'i.4ft).T(4n,4ii)#TH(40,4H),TR{4n,AO).TI(4'.A|i) 

nOMMON/C4/CLFN.NGSXRN,NPY,  SOUND , NMACH, NFREO, MAUG, H I OmCX ,  RHO» 
ft  NMOOFS. LCOLL,LPRWSH,LPRCO»  1 1 Y  *  IIX.NSDRP*  ISftLAT  .  f  W.FC, 

/  Nr.'OI  S,NOM!T.MACH.XCO!  L.YCOLL»PI*U,QWCX,CXHm,  IMOi).  i ROM» 

H  EM.FK#B/»NH!X*NCIX.CHON>NWCY.  IFR.Ei  » E2»  QWY .  QWMX* 

9  SN.WflO,NIY.NWCX/NCnX.NWPX»NCPX,KXWfNG.NYWlMG»NXCS»NYf*S 

COMMON'' Ob/LPHNCH.KF 

EQMIVALFNCF  (AA,AIC)»  (  A  ,  UASH ) ,  <  FH ,  T  ) 

c  » » » 

r  *  *  * 

WRnF(A„6M 
A  F ORMA  T  (  1  Ml  ) 

1  CAI  I  KFIlA 

F W  =  P*NW  I  X  ‘  ’ 

FC  *  P*NCIX  *  1 
OWWX  =  p.il«Pf/FH 
OMC  X  =  P.'l*pf/FG 
QHY  =  PI/FLOAT(.-»#N!Y) 

CAI I  GFO« 

CAI  I  FRAHP(NI  Y.NMr:X,NCCX#NXWING,RYHlNG*NXt:S,NYCS,?*HfJO.SN) 
NRS=N|Y*(NWCX*NCCX> 

NCS=NXWING*NYW1NG*NXCS»HYC5 
no  *  n-j  nit=i  .nrs 
no  phi  hi r=i . ncs 
?!>('  TFHP<NIT.HITI  =  TR(MIT.MIT) 

CHI  TRAMP(NIY  . NWGX , NCCX . NXW I  NO , N YM I NG » NXCS , N YCS , 1 . WpO . SN ) 

no  /in  n i  i  =  j  .  nrs 

00  .01  mi  T=1  .NCS 
TR(Nl T.MlT )=TFMP(MIT.MIT) 
no  111  it  HACH  =  '.NHACH 
MACH  =  MACH 
FM  :  FMACH(MACH) 

0  =  F  M»  SOHNIi 
ROII  -  WOO/II 
RP  I . H  -  I M»FM 
CAI  I  K 0 II T  (  I  ) 

II  (ICOIL  .NF.  0)  CALL  KOUT(  '-> ) 
no  .  llii  |  F  R=  I  .  Nf  RFO 

I  f  R  =  I F  R 

II  <  K  F  .FQ.  i)  t  RFQ(  It  R)  =  FREQ<  It  R)«U/(HHO*?.'i»Pi  ) 

FK  =  .  .  il  »P  I  »FRFO(  IFR)*MHO/IJ 

F  NSURI  =  1  (WIND)  OR  P  (CONTROL  SURFACE) 

NCX  =  NWC  X 
MOM |1  -  1 

N  MO  Hf-  S  -  NXW  I  NO»NYMI  NO  ♦  NXCS*NY(?S 
HAIM.  -•  NCOI.S  »  NHODFS 
00  -  .1-1  .  NCOI  S 
II  (.1)  -  n 
DO  '•  K  -  I  .  MAOC 

A  AA(.I.K)  =  C7FR0  Q 

I  ROW  =  I  -°9- 

nn  i b  NSURF  = I .  ' 


MSHRF  =  NSURF 
HO  14  ! Y  = 1 #  N I Y 
IIY  =  ! Y 

IF fM)HIT-MQMfT,LT.fl)  00  TO  > 
Iff!  Y-NOM{HOH  I  T  )  .  FQ.  n  )  GO  TO  ' 

7  YC-OLL  =  SH»Y(  1 Y  > 

no  i?  lx=t,Nrx 

HX  =  IX 

xcni  L  =  XSd.NSlIRF.  IX,  IY) 
cai  i  conn 
no  ‘.fi  M=  1  »  NHOlif  S 
SR  --  I  R(  I  ROW.  Hi 
S!  -  II  ( fR0«.K1*FK/HR0 
HNf  =  NCOLS  *  H 
vi>  WASH  ( MNC )  =  CMPLXISR.SI) 
no  f.n  n=i  , NCOf  s 
fcf  HA  I  I  CGRFI)(A.N.N> 
f  ROW  =  IROW  »  » 

1?  OONTINUF 

00  TO  M 

iS  MOM  IT  =  MOM  IT  *1 
1 4  CONTINUF 
NCX  =  NCCX 
lr>  CONTI  NUB 
CAU  XL  SO 
CM  l  AICS 
CAI  l  KOUT(fi) 

IF  (l  PUNCH  .  N£  .  'll  CALL  KOIJT  C  /  > 
t  -10  CONTINUF 
00  10  1 
FNO 
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CKfliA 


I  1 

i  '> 


Kf  DA 

SUBROUTINE  KTDA 

COMPLEX  A,  AA, ANH.CZERO.WASHiAlC 
f)  J  Mb  NS  I  ON  A  I  C  {  *  i) »  4ft  ) »  MASH  (  40  > 

COMHON/Ol /A<«"  »,  AA(  40  .  «f) ) ,  ANM<  <0 . 40  > »  CZERO 

OOHPON/CP/HKERf  5>0  )  *  ZKER(?U  >»FMACM(M.FRPG<1  »  >,NQM<5  )  .  JL  (*0  ), 

I  H00R(6> .ZC0R16).  WXOHH<  t  l),WBCN<  jl  )  >  WB I  N(  1  jl  >  •  MT  (  VA  )  , 

V  XM‘.) ,  YE(.l>,UX(  1  0).UY(  1  1  >  ,51X4  41,  .  x  )  ,SC*<  *  «  .  c  » * 

a  El  a  <  •  n 

COHMON/n.l/Y(  •i).XAfC(in#1u,21,VAIC(llt»R)«B(4ii»40)#R(«U»4fi)> 

1  n(iii.4fl)*T(4n,4il)«TM(40,4il),TR(40«4(i),TI(4q,4u) 

noMMON/d/CLFN.NGSKRN,  HP Y.  SOUND,  NMACH,  NFREO,  HAU8,N I OnCX,RHO, 

6  NHnilFS.I.COLL.LPRWSH.I  PRCO,  1 1  Y*  !  IX,NSURF,  ISftLAl ,  FM»FC, 

/  NOOI  5»NijMIT,MACH*XCQl.L«YnGLL*P!»U»GWCX*CXMN>  IHOD.  1RQM» 

»  EM,fcK#R/*NW!X*NClX.CRON#NHCY» IFR.K1  * F2 , OWY . QVfWX, 

SN. WBO.NI Y.NWCX,NCnX.NWPX»NCPX,HXWlNG,MYWlW6»HXCS»KVnS 
rOMHON/n^/LPHNrH.KF 
Ron  I  valence  <aa,aic>. <a,wash> 

FOPMATC*Fl?.H ) 

FOPMKMl?) 

RE  A  !!(•».  H  )(XFC  1).  1  =  1  .*>) 

REAIi(-,,i1  >  (  YE  (  l  ) ,  !  =1  SOUND 
RHOri  .n 

RFAI.  (S,^)  NMACH.KF  . NFREO « LCOLL . LPUNCH 

REAP  (S,K)  NWrX.NCnX.NlONCX.N!  Y.ISOLAT 

INWlSnl 

NWPX=NWCX 

NXWlNO=NWCX 

ncpx=noox 

Nxns=Ncnx 


NWO  Y  =  N I Y 
NPY  =  N|  Y 
NYW1 NO=Nl Y 
NYC<;  =  NJ  Y 

NCOl  S=NPY»<NWPX+NOPX  ) 

noix=nccx«nionox 

NWl X=NWCX»N1 ONCX 

WWTS  =  NWCY«NWCX  *NIY*NCCX 

no  «n  j=i ,nwts 

All  u  1  (  1  )  =  i .  n 

IK  IMWTS.NE.'I  )  R£AD(‘>.11)  ( WT( l > , 1=1 ,  NWTS  ) 

NCO|  S  r  NPY  »  ( NWPX  ♦  NCPX) 

NOIX  =  NCCX«NIONCX 
NW  1  X  -  NWCX*NIONOX 
NOMlT  =  H 

no  *  h  i  >  *i 

A  MOM  <  I  )  r  II 

II  (NWCY.Of  .  N  f  Y  >  00  TO  *, 

NOHIT  =  NIY-NWPY 
PI  A  l‘  (  •>  .  IV)  (  MOM  <  I  I  »  |  =  I  .  NOMI  T  1 
'•  WIADC..II)  ( f' MAORI  I  ) »  |=1  »  NHACH  ) 

00  /  l-l.NMAOH 

II  (I  MAORI  I  ) .  I  E  .  II .  1 9  >  00  TO  / 

W  R  I  If  {t  ,  1  .1 ) 

K  FOPHAM/IM  A  HAOH  NUMBER  OREATER  THAN  n.yu  HAS  BbrN  RFAI)  IN - 

ICASI  fFRMINATin) 

OAII  f  XI  T 
/  OON  I  I  NIIF 

PI  A  n  <‘,11)  <fPEO( I). 1=1 .NFREO) 

PI  A  (.(•>,  i  1  )  (  YA  JC(  I  .  i  ) .  1  =  1  ,  NYWING) 

RE  A  It  ( *i .  I  I  )  (  YA  1  C  (  I  .?).  UI^NYOS) 


IBMWtfiWW  -  -  - 


'SJ4*  ,W0,'J=1  .WYHIKGJ 

YA^ls^i?  ; ;i:;  •"««“ 


***  PH  =  FMACH<1 ) 
^Al  I  K0UT(1) 
SAU  KO»>m> 
RETURN 

eno 


00  70  06 


CKOUT  Kour 

SUBPOUT ! NE  KOIITUND) 

COMPLEX  A,AA.ANM.CZERO»WASK»A!C 
0 1 HF  NS !  ON  SURF  (2,2),  XPR( 50  ) , A I C ( 4  0 , «0  > »  WASH ( *  0  > 
COMMON/Cl/A(8il),AA(40,80),  ANH(40.4f>),CZERO 

COMMON/C?  /  HKFRC20  ).  ZKgR(?t>  ).FHACH(  A)  .FREOdO  >>NOH  (5  > .  1L(50  )» 

1  H00R(6)»ZC0B(6),WXCMM{11>>  WBCN( 11  > , WB IN< 11 ) . HT < 90  ) , 

2  XE(5),YE{3)»UX(lfl)»UY(tO)>WXINMdi)>SIX(40.2)»SCXdO»2), 

3  ETA ( 1 1 ) 

C0HM0N/C3/Y< 1 l),XAlC(in,ln»2)»YAIC(l0,2),B(40,4Q)»R(40»40)» 

1  C(4i|,40)»T(40»40)»TM  (40,40), TR(40»4fl),TIMn*40) 

C0HM0N/C4/CLEN,NGSKRN,NPY»S0UND»  NHACH,  NFREQ,  MAU6,  NI  0»)CX,  RWG* 

6  NMOI)ES,LCOLL.LPRWSH,LPRCO»  IIY,  !  iX.NSURF >  ISnLAT, FW» FC . 

/  NCOI  S»NGMlT»MACH»XCQLL.YCOLL»Pl.U,QWCX,CXHw, I  MOD. I  ROM, 

8  EM.FK,B2,NWIX,NCIX,CB0N,NWCY, 1 FR, El , E2, GWY . OWWX, 

9  SN.WBO»NlY,NWCX»NCCX.NWPX*NCPX,NXWlNG»NYWIwO*?iXCS>NYCS 
C0MM0N/C5/L PUNCH, KF 

SOU  I  VALENCE  ( A A , A fC ) , ( A . HASH ) , < XPR . I L  ) 

BATA  (  SURF  (1  .  I  ),  I=1»?)/8HWING  ,8H,TAfl  / 

BO  TO  (10.2H.  If),  40, 50,80, /0>.  JND 
C 

1 n  XV=XE(5  )-XE( 4  ) 

XX  =  XE  (3  )-XE(  ">  ) 

AW  =  2.u«(XE(3>*YE(.>)-0.5»XE(2)*<YE<3)-YE<2))  > 

AT  =  2.ll#XV*YE(3) 

WRITE (A, 1 l ) EH , SOUND, RHO, XE ( 1 > » XE ( 4 ) #  XE( 3 ) , X V , YE ( 2  >  #  YF ( 3  ) , YE( 3  > » 

1  YE (3), XX, XV, AW, AT.NWCY.NI Y . NWCX , NCCX, NW I  X . NC l X » NPY , NPY » NWPX# NGPX 

11  FORMAT  ( 1  HI  /////  .VX.41HHUGHES  AIRCRAFT  CO.  SUBSONIC  AfC  PROGRAM 

1  ///3/X.30HFI IGHT  CONDITIONS  AND  GFOMETRY/J H0//15X,  i 3HHACH  NUMBER 
?  =.F8.5,4X,15HSPEE0  OF  SOUND  =F10.3>4H  L/T .  4X ,  4HRH0=  ,  El  4 . 8//1H0/ 
3<MX,4HW!NG.  )  0 X . 

3  4HTAIL///22X.16HL.E.  STATION  (L > , 2F22.3//22X, 1 6HR00T  CHORD  (L), 

4  2F22.3//  22X.16HL.F.  SPAN  ( L > , 2F22 . 3/ /2 2X , 1 6HT . E .  SPAN  ( L  )  , 

5  2F22.3//  22X.16HTIP  CHORD  ( L  ) ,  2F22 . 3//22X ,  J  6HT0T  *L  AREA  <L*L) 

6  2F22..1//22X.  l^HSPAN  COLL-  STA.  .  II  9,  122, //2?X»  16HCHnRD  COLL.  STA. 
/  II 9,  I22//22X.16HCHORD  INTO.  STA . ,  1 19 , 1 22/ /22 X . 16HSP»N  PRES  MODES 
8,  Ii«»  l??//22X, 14HCH0RD  PRES  MODES , 1 19 , 1 22 ) 

If  (f  MACH(MACH).LE.0.9»>)  GO  TO  1  5 
WR 1 TE ( 6 , 1 4  ) 

14  F  ORMA  T ( 92H  A  MACH  NUMBER  GREATER  THAN  0.95  HAS  BEfcu  USEO - 

1 USF  CAUTION  IN  APPLYING  CASE  RFSULTS ) 

15  I  F ( NOM I T  *  EO . 0 )  RETURN 

WR1 IF (A. 1 2 ) <NOM( I ) . 1=1 , NOM IT ) 

12  FORMAK1HII.  I'*X,51HTHE  SPANW1SE  COLLOCATION  STATION(Sl  OMITTED  ON  W 
I  1  NR  .  9  1 5  ) 

RETURN 


20  NCX=NWCX 
N I  X  =  NW I  X 
no  15(1  NS  =  1  .> 

WR1 IF (6. 22) ( SURF ( I , NS ) , 1=1.2) 

22  FORMAT ( 1  H 1  , / / / /30 X , 50HHUOHES  AIRCRAFT  CO,  SUBSONIC  A[C  PROGRAM  (CO 
INT-ll)  ////// 

22HX.53HUNSTI A 0 Y  AERO  COLLOCATION  STATION  COORDINATES  ON  THE  2A6/1H 
3 0 , 1 2 H  S  STA  N0.7X.2HYC»8X,/X,11HXC  VALUES--) 

no  121  I  Y  =  1  ,  N I  Y 

YC=WBO*SN*Y( 1 Y > 

no  1211  JX  =  1  ,  NCX 

120  XPPf IX)=WB0*XS( 1 ,NS. IX, IY) 

12.1  WRITERS,  124)  IY,YC,(XPR(IX).!X=1.NCX) 

-93- 


124  FORMATOHU,  I1v,5El7.6/(lH  .  29X .  4E17 . 6  > ) 

WRITF(6,105)  (SURF( I, NS). 1*1.2) 

1  (1*5  FORMA T( 1  HO, 2 4 X.39H INTEGRATION  STATION  COORDINATES  ON  THE  2A6/1HU, 
ll?H  S  STA  N0,7X»2HYI,8X,7X»1lHXI  VALUES--) 

00  106  I Y  =  1 »  N I Y 
YI =W80*SH»ET  A ( I Y ) 

00  1 ?6  IX  =  1 .NIX 

176  XPP( IX)=WB0«XS(2.NS, IX, JY) 

106  HRITE(6,124)  I Y, YI * ( XPR( !X>, I X=S .NIX) 

NCX=NCCX 
N l X  =  NC I  X 
1*0  CONTINUE 
NXS=NXWIN0 

nys=nyhing 
no  ?no  ns=i ,7 

WRITE  (6,201  )  ( SURF ( I , N$) » l *1 »  ? ) 

?Ul  FORMA T ( 1  HI , / ///3flX* 58HHU8HES  AIRCRAFT  CO.  SUBSONIC  A?C  PROGRAM  (CO 
1  NT-0 )  ////// 

2?8X,  4.SHA !  C  COLLOCATION  STATION  COORDINATES  ON  THE  2A6/1H0, 

319X,  4HYAIC,13X,13HXAIC  VALUES--) 

DO  ?n2  IY=1,NYS 
YC=YAIC( IY,NS) 

20?  HR  I TE ( 6 , 203 )  YC, ( XA IC( I X, I Y, NS ) , I X=l, NXS ) 

NYS=NYCS 
200  NXS=NXCS 

?03  FORMAT(1HO,17X.5E17.6/(1H  , 29X , 4E1? .6  ) ) 

RETURN 

C 

10  no  34  NS  =  1.2 

WR!TF(6,?1)FREQ( I FR ) . NMODES, EK , EM 

21  F0RMAT(1H1,31X.4?HMISSILE  SUBSONIC  AIRLOADS  PROGRAM  fC0NT-D)//lH  / 

1  9X »  27H0SC I  LI.  A  TORY  FREOUENCY  ( CPS  ).  FI  2 .5, 14X ,  I  ?,  1 5H  rOLL.  STATIONS 

2  /I  HO , 8X» 3UHRFDUCED  FREOUENCY  ( SEMI-CHORD ) , F9. b, 14X . 23HFREE  STREA 
3H  MACH  NUMBER,  F9.  ,1/3  M  ) 

WRI TE (6, 31 )  IMOD 

31  FORMAT(31X,40HPRESSURE  COEFFICIENTS  FOR  COLL.  STA.  NO. I3//19X, 1HI 1 
21X.10HR  COEFF ( I )1 2X, 10HI  COEFF(I)  9X,9HSPAN  MODE  3X,lOHCHORD  MODI) 
HRITE(6,32)( SURF (KI,NS)»KI=1»2) 

.5?  F0RMAT(1H0,9X,?A6//) 

GO  T0(?,3),NS 
?  NL  =  NHPX 
ML  r  NPY 
I K  =  l 
GO  TO  4 

3  NL  =  NCPX 
ML  =  NPY 

IK  =  NWPX*NPY*1 

4  no  6  I MM= 1 , ML 

no  6  INN-1  . Nl 

HR  I TF(6,33)  IK, ANM( IK, IMOD). I  MM. INN 
13  FORMA T(1 HO. 1 1 9 . 1P7EP2.S, 2  1 1 3 ) 

6  IK  s  |K  ♦  1 
34  CONTINUE 
RETURN 

0  ••«•*»•««•••••»•«•••»••••• 

40  HRITF(6,4l> 

41  FORMATMH0.20X.  58HFRROR  IN  INPUT  DATA  (NO  TAIL)  REOU IRES// » 21X,  19H 
1  TERMINATION  OF  CASE) 

CAM  EXIT 


•>1  FQRmAT(1Hi)»/0X»  63HNUMBER  Of  COLLOCATION  OR  I NTFGRAT  ION  STATIONS  OK 

1  PRESSURE  TLRM5///1  X»2*>HEXCEEDS  ALLOWABLE  MAXIMUM/ / /  35X# 18HCASE  IS 

2  TERMINATED) 

CALL  EXIT 

0  #»»t»##»##i«f#«#*»»«»t####»#»#*»*«»«****»***»*»**f*****4***#***,*# 

60  VEL=FM*SOUND 

0=fl.5«RH0»VEL#*2 
EKt  =1 . fl/EK 

REFC=(XE(3)-:<E(t))/2.0 

WRITE  (6,220)  FREQ< IFR ) , R£FC, EK , EK1 , EM, U, RKO. Q 
??0  FORMAT  ( 1  HI ,  31  X, *»nH  HUGHES  AIRCRAFT  CO.  SUBSONIC  AIC  PRO0RAH  (CONT 
1-D)////9X» 28H  OSCILLATORY  FREQUENCY  (CMS ) .  4X»  1PE12.5  . /I  HO  , OX ,  15HRE 
2FERFNCE  CHORD.,  4 X,  IPE12. 5, /I HO , 9X  . 30HREDUCED  FREQUENCY  (REF.  CHORD) 
3.4X.1PEl2.5,/1Hfi.9X,29HRE0UCED  VELOCITY  (REF.  CHORD ).4X,1PE12.5, 
4/1  HI)  ,9X,23HFR£F  STREAM  MACH  NUMBFR ,  4X,  1PE12  .5 ,  /  J  HO ,  9X,  20HF  REE  STRE 
5  AM  vELOClTY,4X*lPE1?.*i,/lHO»9X»'/HDFNSITY,4X,OPFb*2»/lHO,9X»33HDYNA 
6MJC  PRESSURE  ( 1 /2*RH0*VEL»»2 ) , 4X, 1PE1 2.5 , //// ) 

WRITE(6»2?1 ) 

2/1  F0RMAT(///36X,  14HAER0D YNAM l C  INf  LUFNCE  COEFFICIENTS# //4X,2HRL#1 OX, 
1 2HIM, inx»  2HRL,  1  OX.2KIM. 10X.2HRI  » t  0 X# ?H I M , 10 X , 2HRL ,  1  0 X » ?H I H, 10 X , 2HR 
?L,H<X,?HIM,  /  ) 

NROWS  =  NYWIN(;*NXWING«NXCS*NYCS 
no  222  NROWsI .NROWS 
WR!T£(6.223)NR0W 

WRITE (6, 224 )  ( AlC(NROW,NCOL).NCOLs1 #NROWS) 

22.1  FORMATS  5HROW  =  12) 

224  FORMAN1P10E12.4) 

222  CONTINUE 
RETURN 

C  #«♦«»##»»»*»»##«•«#•«»♦»»♦#«#»««#««*•»•*•»«*•*•»#»»•»#•»**•**••**» 

/0  NW=NXW!N6»NYWING 

nc=nxcs»nycs 

NT=NW»NC 
NW1 =  NW  +  1 

GO  TO  ( 81 »  82 • 83 , 84 ) »  LPUNCH 
81  CONTINUE 

DO  .101  1  =  1, NW 

PUNCH  85,  <  A  I C ( I , J ) , J=1 , NW ) 

301  CONTINUE 

88  FORMAT  (1P6E12.8) 

RETURN 
8?  CONTI NUF 

DO  .102  I  =NW1  .  NT 

PUNCH  85.  ( A  I C ( I ,J),J  =  NW1 , NT ) 

302  CONTINUE 
RETURN 

83  CONTINUE 

DO  i03  |  =  t , NW 

PUNCH  85.  (  A  I  C  (  I  , .)) .  J  =  1  .  NW  ) 

3iM  CONTINUE 

DU  304  | =NW I . NT 

PUNCH  88.  ( A  I C (  I , J  ) , J  =  NW1  , NT ) 

304  CONI  I  NIIF 
RETURN 

84  CONTINUE 
DO  308  f  =  1  .NT 

PUNCH  85,  ( A  I C ( I.J), J  =  1 .NT) 

106  CONTINUE 
1ft"0  RETURN 
FND 


-95- 


tuics 


SURKOUYSKfc  AfCS 

OOMPLFX  A»AA.ANH,nZFRO,WASH»Ain 
0  IMF  NS  JON  A|C(4ft.un).HASH<4U)»FH<4fi#4U) 

nOMHON'CI  /  A(Ri<  ),  AAHn,A0)»  AMH(*u.4n  ),CZER0 

COMMON 'C?/KKFR< 20  > ,  ZKER(  ?fl  >  >  FHACH<  *  ) ,  FRfcQt  1 1I  ).N0H(5>)  .  1L(*U), 

1  HC0R(6) .ZC0R(6)»MXCMN(1 l)»WBCN(11  ) .  WB I N( 1 1 > , WT <  V 0  ) , 

?  XF(r})*Yf;(.0*IIX(lfl)>UY(iU)»UXlHN(U)*SlX(4b  . 2 > . SCX ( 1 0 . 2 > , 

3  E  T  A  ( 1  n 

COMmON/C3/Y<  H),XAlC<in,ll!.2),YA!C<in,?).B<4  0.40>»R('ili»40), 

1  C(4ii.4n),T(4n,4M)<TH{4n,40),TR(4n,40),TIC4o,4|i) 

C0MH0N/C4/CI  FN.  NOSKRM.WPY,  SOUND,  NMACH.  MFRhO,  MAUG.  N! 0>jCX ,  RHO. 

'»  NMOOFS,lCOLL»LPRWSH»LPRCO»nY»llX»NSURF»  I  SnLAT  . F  W»  FC . 

/  ncols, nghi t « mach. xcoll  .  ycoll*  pi  »  ii,  qwcx.  cxhn, i mod,  i row. 

B  FM.FK,R^,KW!X,NCIX.CHON,NWCY, IFR.fc’ ,  F2.UWY . UHWX, 

V  SN. WRO.NIY.NWCX.KCCX.NWPX.NCPX.WXHfWO.NYWlwG.NXCS.NYCS 

FUUI  VALENCE  UA,AIC>.(A,MASH>»(FH.T> 

WCOl  S*NPYMNWPX*NCPX  J 

HROWS  =  NXWIHG*HYW!RG-»MYCS«NVrCS 

CALI  FORCE 

00  /O'!  I  =  L  »  NROWS 

no  ''Oil  J=1  ,  WOOL S 

SR=FM< 1 , J) 

SI  = 

?»n  AAM  ,  J)=CMPLX(SR,  SI  ) 

00  -.0.1  1*1.  NROWS 
no  *11  J=l,  NROWS 
A(j)=(n.n,n.n> 
no  s ii  k=i  ,Ncm  s 

?|0  A( J)SA( J)-AA( | .K)*ANM(K, J) 
no  •/?•>  I  * L  ,  BROWS 

A1C(!,I.  )  =  <YH  -.)-Y6n  >  >»A<L>/(/.M«EK»«2> 

YtlO  CONTINUE 
RETURN 
ENO 
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CCGRFD  CORED 

SUBROUTINE  CORF  0 ( V,  JR* !C  ) 

COMPLEX  ANH#  CZFR0. A  1C . WASH 
0 1  MENS  I  ON  A  I  C(  4|) ,  Hfl  ) .  WASH(  40  )  #  V(  ?»  1  ) 

C0MM0N/C1  /A(?,«l)).AA(2.40,fiO>, ANH<40.40>.CZERO 
COMMON/C?/HKFR(2n).ZK6R(/?U)*FMACH(A)»FREQ(10)»  N0H(5  ) » IL (So ) . 

1  HC0R(6)»ZC0RI6)»HXCHN(11)»WBCN(11  > .  MB  !NU  3  ) ,  WT<  90  > , 

2  XE(0,YE(3),UX(l<n«UY(1u>>WXIHM<1  l ) ,  SIX<4«  .  2  ) ,  SCXt  1  II .2  ) » 

3  ETA ( 1 1 ) 

C0HM0N/C3/ Y<  1  l)»XAJC(1A,t(j.2)#YAfC<10»2)»B(4r'»40)»R(40*40), 

1  C(4H.4O)»T(4O»4H)»TH(4O#A0),TR(4fl*40),TI(4'»40) 

COMMON/C4/CIFN.NGSKRN»NPY#SOUND»NHACH.NFR£Q»MAUG»NIOnCX,  RHO» 
ft  NMODFS,LCOLL»LPRWSH.LPRCO>  I IY,  I  IX.NSURF* ISftLAT.FW.FC- 

/  NCOLS.  NOHIT, MACH.  XCOLL .YCOLL. P|, U.OWCX.CXMn.  I NOD.  IRON, 

tt  FM.FK.R'.NUIX.NCIX.CRON.NWCY.  IFR .  EJ  #  E2#  QWY .  QMnX , 

9  SN.  WRO.  Nl  Y.NWCX.NCCX.  NWPX. NCPX , NXW l NO, NT W I  «*G  •  NXCS»  N  fCS 

FOIII  VAI  f  NCR  (AA*AIC).(AtWASH) 

RMN  =  SORI  (AAO  »  IR,  IC)*«?  ♦  AA  (  2  .  I  R#  I C  )  •*?.  ) 

IF  (  A  A  ( ? »  IR.  IO.LF.E2)  GO  TO  3'i 
CR  -  A  A  ( 1  #  IR.  IO/RMN 
Cl  r  AA<?,  IR,  IO/RMN 
DO  .  A  N- I C#  M A UO 

T I  =  CR«AA(1.IR.N)  ♦  C I  * A  A ( 2 . 1R.N) 

A  A  (  /  »  IR, N  )  =  C  R  *  A  A  ( >  » IP.N)  -  C I »  A  A  <  1 ,  I R,  N  > 
vn  A  A  (  I  »  I R ,  N )  =  T! 

>n  RAN  =  SORT! V< I . IO*«?  ♦  V(?»!C>»*2>- 
IF (RAN.LF.E2  )  GO  TO  ftA 
RAN  =  SQRTtRAN**?  ♦  RHN**2> 

CR  =  V < 1 . !C  )  /RAN 
O  =  V(2.IC)/RAN 
RMN  =  RMN/RAN 
DO  ‘  «  N= I C »  MAUO 

AIR  =  RMN*AA( i . IR.N)  ♦  CR*V(1,N)  *  CI»V(?.N) 

All  =  RMN*AAC/,  IR.N)  ♦  CR*V(?,N>  -  Cf«V(1»N) 

VR  -  RMN* V ( 1  . N )  -  CR*AA(1 , IR.N)  ♦  C I  •  AA '  ? ,  I  R . N > 

VI  t  RMN*  V  <  2  .  N  '  -  CR*AA<?»  IR.N)  -  C I  * AA ( 1 , I R . N > 

A  A  < l  .  IR.N)  =  AIR 
hue/,  IR.N)  =  All 
V  < 1 , N  )  =  VR 
•>n  V  ( ?  ,  N  )  =  VI 
ftn  RETURN 
FNP 


CFOHCF 


StJRROuTINE  FORCE 
n  IMF  NS  ION  FK(4fl.4«>> 

COMhON/C3/Y( 1  i  >»XA?r(l»»t«.2)#YAlC{lfl»?).8{4J»/4U)»R(  »lt.  «}()>, 

1  C(4ii.4n).T{4n#4il),TH(4H.40),TR(4n,4|l),Tl(4n,4u) 

COMMON /C2/HKFR(2n ),2KER(?U > . FHACHC * > » FRfcQI 1 0) ,NOH(5> . IL(5U ). 

1  HCOR 1 6 ) »  ZC0RF6 ) »  WXGMNt 11).  R3CN( 1 1 ) . MB ! N(  1 1) , MT  ( 90 ) . 

H  XF(b),  YE(3),UX<  lO>»OYCtfl  J.MX!HI|(11)#SIXC4U.2>»SCX<1U.^)» 

3  ETA(U) 

C0MH0N'C4/CLFN.NGSKRN,NPY.  SOUND#  NMAGH»NFR£Q,HAUG»NIO»iCX,  RHO, 

A  NMODFS.  LCOLL  »  LPRWSH>L.PRCO*  ISY»  !  IX.NSURF#  IS'iLAT.fW.FC. 

7  NCOLS# NOHIT,  MACH,XGOLL»YCOLL»Pl,U»QwCX»CXHw,  I  MOO.  I  ROW. 

8  EM,EK,B2.NWIX,NCJX.CR0N,NWCV.  IFR.E1  ,F2,0MY.0MWX, 

9  SN«W80.NIY»NHCX#NCCX.NWPX/NCPX*NXMlN6*NYWl*i6»wXCS»NYnS 
FO»H  VALENCE  <FM,T> 

MROwS=NIYM NUCX  +  MCCX ) 

MCOI  S=NPY»(NMPX»NCPX) 

00  1  On  I  =  1  »  HROW5 

no  !«»  j=i. Knots 

t  '1ft  FM(  J  .  J)  =  '*.n 

n  *»*  BEG I N  TO  ASSFMRLE  FHUROW,!COL>  MATRIX  STARTING  WITH  WING 
IROW=l 

no  .'OH  1  =  1,  MY 

IF  {  I  .EQ.  1  1  GO  TO  !«•» 

IF  ( i  .EO.  N I Y  >  GO  TO  110 

SNt  -  <  •!.•»•<  YA ICf  l- 1#  1  )*YAIC<  I  •  1  >  >-YF(1  ))  /  (YF(  M-YEC1  >1 
SNU=(il.*;*(YMC(  I,  i  )*YAlC(I*l.i>).-YF(1  >)/{YF(i)-YE(1  M 
GO  TO  IIS 

I  lib  SNI  rP.n 

SNU=  ( il . 5M  YA  ICC  l » 1  )*YAIC(^»l))-YE(U)/(YE(3)-YeU  )  ) 

GO  70  its 

II  n  SNI  =  (Ii.(j*(YAICINIY-»  ,1  >  ♦  YA  I C  C  N I  Y  -  1  )  I  -  YE  ( 1  )  )/( YF(3  >-Ye{  1) ) 

SNIM .n 
1!8  CONTINIJF 

no  /no  j=i,nwcx 

IF  ( J  . EO.  1  >  GO  TO  1  ?n 
IF  (J  . EO.  NHCX)  GO  TO  l?* 

CNI  =  (  •» . S *  <  X A  I C; <  J.  I  .  1  MXAICC  J-l  •  1  •  1 )  >-XE<1  )-«.!>*(  XE ( 3  > -XE  ( 1  >  >  )/ 
lIll.l>*<XF(.i>-XE(t))> 

CNII- fl.bMXA  I  C(  J*  i.l.l)  +  XAlC(J.).l))-XE(1)-0.b*(XE(J)-XE(l)))/ 

1  (U.‘.#(XE<3)-XF(1 ) ) ) 

GO  TO  130 
1  'f  CNI  =  - 1  .11 

CNII-OI.SMXA  ICII  .  I  .  1  )  *  XAICC2.  I  •  1  )  )-X6(1  1-0.5*  (XE(3)-XE<  J  ) ) )/ 

1  (lt.H»(XE(.f)-XF(l  >)> 
no  TO  130 

CNI  -  ( •).*»•(  XA  1C( NMCX-1  .1,1  MXAincNWCX,  1 .  J  ) ) -XF  ( 1  )- . *# f  XF < 3  ) -XE ( 1  )  ) ) 
I /<0.*>*<XF<3>-XF(t  )>) 

CNHr  < | .0 
1  in  CONI  INIIE 

I  not  -i 

no  ;-nn  k  =  i  .npy 

CAII  M l N  TS( K • SNL . SNU# FS ) 

no  Mill  l.  =  l  *  NMPX 

CAM  MINTCd  .CNL. GNU, FC) 

FH< IRON. ICOL ) - F  S • F  C 
5Mi  0  I COI  =  I  COL  ♦  1 
Ion  I ROw=  f  ROW ♦ I 

C  *•«  ASSFMBIF  CONTROL  SURFACE  CONTRIBUTION 
00  /.fin  I  =  1 » N I  Y 

IF  (I  ,FO.  I  )  GO  TO  *» it *»  -98- 


IF  {  I  .EO.  NIY)  00  TO  Sifl 

SNL  =  (lI.b»(YAIC<  I-L»?)«YAIC<1  .  <>>)-YF<1  ) )/ (  YFO  >-Y£<  1  >  > 
SNU=H).S»<YAIC(  I.  7)«YAIC(I  ♦!»?>  )-YF(1  )  )/  ( YF(  >-  YE  C 1  >  t 
00  TO  S?S 
SuS  SNI  =0.0 

SNU=  (O.S*(YAIC<l»^MYAICI?*2))~YF(1  ) )/  ( YE(  3  )-  YE  ( 5  ) ) 

00  TO  StS 

SJ  ft  SNLMH.h«(YA1C(N!  Y-1  » 2  HYA  IC(  H  J  Y.  ?. )  )-YE<  1  )  )/<  YE  <3  >«YP<  1  )  > 

SMU=1 .B 
SIS  CONTINUE 

00  MHJ  J=1 »  NCCX 

IF  ( J  .EO.  1)  00  TO  7?0 

ir  (J  .EO.  NCCX )  GO  TO  7?5 

ONI  =  (i|.*>*(XA|C(  J.  »  .?>+XAfC(  J-l.l  .2)  )-XE<  4  J-il  .5*(XE<!>  >-Xfc<4  > ) )/ 

1  («.'.*{  XF<S)-XFH>>) 

CN»=(».S*(XA  ?CF J* J  » J ,?)*XAIC( J. J .?) )-XE(4)-0.S*(XeC&)-Xfe(4)  >  )/ 

1  C  «.».*<  XF(*)-XF(4>>> 

00  TO  /.ill 
T/W  CNI  =-t  .1* 

CNU-(it.S*(XA!CIl »  I  >7  )  *XA  JC(2  •  i  •'/-)  >-XE { 4 )-o . { X£< b >-xE ( 4  >  >  >/ 

1  (i*  .S*(XE(S)-XE(4)  )  ) 

00  TO  7.5  0 

7^S  CNL=(fl.b*(XAIC(NCCX-1  .1  ,  ?  )♦  X  A  I C  (  NCCX.  1 . 2  )  )-  XE  (  A  )~0  .  b#  (  Xfc  ( t>  )-XE  (  5  >  i 
J. ) /  ( i< . S*  (  X£< b  )-XE(  4  > ) ) 

CNU=1 .0 
7.5  0  CONTINUE 

I  COl =NPY*NWPX  *  1 

no  »  flu  K  =  I »  NPY 

CAM  MSNTS(K.SNL.SNU.FS) 

no  ^ mi  l=i,ncpx 

CAI 1  MINTCIL.CNL.CNU.FC) 

FM  (  IRON.  I  COl.  )=FS»FC 
non  I  COl  -  I  COl  ♦  1 
6im  I ROW= I ROMM 
RETURN 
FND 
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coord  cord 

SURROUT  I NE  CORD 

COMPLEX  A,AA,ANM.CZ£RO»WASH»AIC,AK.H?,TRM,Rl..CRNL 
DIMENSION  AlC(4i),yn),MASH(40) 

COHMON/Cl  /  A  {  Of) ).  AA  <  40,80  ),ANMHtl.  4ft),  CZERO 

COMHON/C?/HKFR(2n)»ZK£R(2ft>»FHACH(M,FREQ(10)>NOH(5)  .  IL(Su>, 

1  HC0R(6)  ,ZC0R(6).«XCHNni)»WBCM(ll  ) ,  MB  IN<  1 1 ) ,  WT(  9fl)  . 

2  XE(*>#YE(3).UX(10)*UY<10)»WXIMN<11  )#SIX(40.?>#SCX(1U.2'). 

3  ETA(li) 

COMMON/CVY<  1  I  )»XAlC(in,l0»2)»YAlC(10»?),B(4ft»4U)»R(  4 1)  #  4 II )  * 

)  C(^«.4f>),T<4n,40),TH<40»40),TR(4n,40)»TI(4n,4«) 

COMMON/C4/CLFN. NGSKRN. NPY, SOUND, NHACH. NFREQ, MAUG, N l OuCX , RHO, 

6  NMODFS. LCOLL.LPRWSH.LPRCG, I  I Y  # I IX.NSURF. I SnL AT .  f  W , FC . 

7  NCOI.S,  NOHl  T,  MACH.  XCOLL  •  YCOLL >  P I , U,  ONCX,  CXMn,  I  MOD.  I  ROW. 

8  FM.FK,B2.NWIX,NCIX,CR0N,NWCY. IFR.E1 ,  E2 »  GW Y . QWWX » 

9  SN •  MRO.  N!  Y ,  HMCX •  NCCX .  NMPX.NCPX ,  NXWf  N8 ,NYWi NG  >  NXCS*  NYf'.S 
EQUIVALENCE  < 4A, AlC), (A, HASH) 

C  THIS  SUBROUTINE  CONSTRUCTS  A  ROW  OF  THE  DOMNMASH  HAToIX 

C  THE  PRESSURE  SERIES  IS  A  PRODUCT  OF  CHEBYSHEV  POLYNOMIALS  IN  THfc 

NEGATIVE  OF  PERCENT  SEMI-CHORD  FROM  THE  MID-CHORD  aND  PERCENT 
SFM I -SPAN  FROM  THE  ROOT. 

DO  6  JCal.NCOLS 
8  A ( JO  )  =  CZERO 
IC1  =  ft 
NIX  =  NWIX 

OMX  =  -QWWX*SN**?/<8.n«PI ) 

NPX  =  NWPX 

C  THF  DO  14  LOOP  COMPUTES  THE  NON-SINGULAR  PORTION  OF  n(N,M> 

C  DUF  TO  BOTH  SURFACES 
DO  14  NSURF  =  1 

IF ( MSURF . NEcNSURF. AND. ISOLAT.NF.n )  GO  TO  13 
DO  1?  I  Y  =  J.  #  N  I  Y 
ETAI  =  SN*ET  A (  I  Y ) 

ET?  =  FTA ( I Y ) **2 

IF(NPY.GT.l)  CALL  CH6B (NPY-1 . ETA ( I Y ) , UY ( 7  )) 

UY ( 1 )  =  1 .0  -ET2 
DO  A  K  =  ?»  NPY 

3  UY ( K )  =  E  T?*UY( 1 ) *UY ( K ) 

DO  10  IX=1 .NIX 

XI  =  XS(?. MSURF, IX. IY) 

X I D  =  XCOLL  -XI 

AK  =  CRNLIEK.XID. YCOLL-ETAl.EM.BP)  ♦  CRNL ( EK . X I D, YCO)  L*ETA I .  EM . B? ) 

IC  =  IC)  vl 
H 7  =  AK#OWX*OWY 

IF(NPX.GT.I)  CALL  CH^B(NPX-1.-SIX( IX, MSURF), UX(2)  ) 

MX( I  )  a  1.0  -S!X< I  X , MSURF ) 

DO  4  K=2»NPX 

4  UX(K)  =  <1.0  ♦SIX! IX, MSURF) )«UX( 1 )*UX(K> 

C  **  ADO  AN  INCREMENT  TO  EACH  ELEMENT  OF  THE  ROW  FOR  (yl.fcTAI)  *• 

DO  in  NV= i »  NPY 
TRM  =  HP  •  UY(NY) 

oo  in  Nxn.NPx 

A ( I C  >  S  A ( I C  )  *  TRM*UX ( NX ) 
in  IC  =  IC+1 

C  *#  IC  EQUALS  NPY#NWPX+1  AT  THE  END  OF  THE  FIRST  PASS  ## 

I  ?  GONTINUF 

II  NIX  a  NCI  X 

QWX  =  -QWCX*SN**?/<«.(|#PI) 

NPX  =  NCPX 
14  Id  a  NPY*NWPX 

I  ci  =  n  -loo- 


n 


NPX  =  NWPX 

XCOtS  =  XS  < 1 » NSURF  *ItX»IlY> 

Y2  =  Y<IIV>**2 

CAU  CHEBINPY-1  .  Y<  II  Y ) .  UY<  ?  >  > 

OY ( 1 >  =  -2. ft 
DO  1*>  K  =?  #  NP Y 

15  UY  ( K  )  =  -2.fl*Y?»UYCK> 

DO  * fl  MSURF  = 1 »  NSURF 

••  THIS  LOOP  ADDS  THE  CONTRIBUTION  OF  THE  SINGULAR  INTEGRAL 
ALONG  THF  I.TNF  FROM  THE  HING  t«F.  TO  THF  COLLOCATrON  POINT 
IFIMSURF. NE. NSURF. AND. ISOLAl.NP. ft)  GO  TO  23 
IHNSURf.LE. MSURF)  00  TO  16 
UPL  JM  =  PI 
GO  TO  10 

16  XT  =  SCX< 1 1  X . NSURF) 

KPI.IH  =  -ATAN(SQRT(1.U-XT»*2)/XT) 

IF (UPLIH.LT.ft.ft )  UPL I MsUPL I M ♦ P 1 

1  F  OWSNG  =  FLOAT(’*NI Y)*UPLIM/b.« 

DO  '•?  N  =  1 ,6 
1C  =  I  Cl ♦ 1 

C  ••  THIS  LOOP  CONSTRUCTS  D(ft.M)  »  M=0 »  L  * . . . . NpX-1 

VINT  =  UPL I M«  ZCOR ( N ) 

C  =  COS ( V  I  NT ) 

CAM  CHEB  (  NPX- 1  *  C.IJXd)) 

IIX <  i  )  =  1  of)  »c 

no  i 9  K  =  ? » NPX 

19  UX(K)  =  (l.t-  -C  )  *UX  ( 1  )  »UX<  K  ) 

ARG  =  FK*  <  XCOLS  -HXCHNUIY)  ♦C*HBCN(  I  I  Y  )  ) 

IF  f MSURF. EG. P 1  ARO=FK#(XC0LS-CXHN*C*C80N) 

Cl  =  COS ( ARG  ) 

Si  =  SORT ( i .«  -Cl ) 

01  r  CMPLXICi ,-S1 >*HCOR(N) 

DO  > ?  NY=1 .NPY 

TRM  =  QWSNG»UY ( NY ) *D1 

DO  f?  NX= 1 »  NPX 

A  (  I C  )  =  A (  I  C  1  ♦ TRM®UX( NX  ) 

2?  ! C  =  I C* 1 
23  ICl  =  NP Y*NWPX 
<\  ft  NPX  =  NCPX 
RETURN 
FNn 


<*  JC 


rrim 


CRHl 

roW’ttX  runi  T|f1«  i.-RNUCK.S.Y.CH.RX) 

CGMHON/C?/HKFR(  M  ) .  ZK£P{?0  1 ,  FNACHi  A  ) .  FRbQ(  1  (’  ) .  NUM(*> )  ,  I  L  { <>o  ) , 

1  Hi  ORifti  .ZfOR<6KWXCKK(  t  I  i.HBCNCH  ) .  HB!N(  11  i .  U1  (  VIU  . 

?  XF  (*- 1  •  Yfe-(  .«>  *tlX(  1*  >  *UY  f  1  0  >  #ttXIHN(  1  1  >,S1X(4U  /  1  *  SCX  ( 1  u  .  I )  . 

*  FIA( i  1  } 

COK*ON'C«/ClFN.NGSKRN,NPY,  SOUNR,  NHACH,  NFPfcG .  MAUli,  N !  U»  UX,PHO. 
f  HHnaFS.LCOLL -LPRHSH.tPRCO* ! ! Y,  I f X.NSUKF.  I S«LA  1 .  I  W,  FC . 

1  NfntS.HOHIY.HACHfXnGU.YnOt  I..PI.U,OWCX,CXHn,  inn£).  IROW. 

EH, PK , B  * . HH i X » NC ! X • CROH , NMCY  • !FR,pl , ti , QHY . QWkX , 
SN.HR0.Nir.NWCX.HCCX.HMPX.MCPX.HXMiN0,HYHiH6.NXUS»NlfrS 

RsAFS(Y) 

R-VaQ  | 

f!Ki  s  f!K*H 
0  i  «( .  K 
04  •».<* 

04 

H  S>  <«  X  •  H  »*R  * 

S  -  SORT  i  S-  )  , 

Hi  iH’*R> 

H*  OK  3  *U1 

no  f  f  =  i.noskrn 

117  *  lit  *ZKFPf  i  i  J 

•I/1’  r  U7«*? 

0  =  IM'e7KFR{ }  l 

r  -  HKFRC1)  SQRTi » ,«*1I2»  >»UZ*U1 
0  '-I  <  *F«FOStR I 
R*=«.  *  « F»S I Nt  0  j 
V  =  1  -  7KFRI l >.* 

F  -  HKFR(!J»  .  *  V  •  PXP< -OK  «  «V  )  <  SQRTU.««V1 
*n  Ol  rl  1  *F 

0  #  -  O’  *  G  * 

t  s  X'S 

If  U  K  .HI  .  ».  •  I  00  TO  >•' 
r  1 4  =  i .  n 

on  in  *>  < 

'*■»  ri  4  r  CKl  »BI  01  fCK*  > 

*  i  oi=i  Ki  *i;4-fi  *-xs*ons(Ux  ) 

G/-rx  i  «r./*  XS* s i n { ok  i 
KK  OK  «  X 

on  nnsmi 

SI  SINiXM 

OMNI  <*NP|  >(  (rt}»ui  *$l  *G?  u:0»R/~St*ni  )/R<*  > 

PI  TURN 
r-Nn 
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R  R  SL 

FUNCTION  BFSUX) 

!F(X.0T.?.0)  00  TO  Sn 

T  =  X,;5.76 
T  =  T.T 

RSI  I  =  'l.6  +  T*<  fi  .P7flvn»j94*  T*  (  0.5  1  ->9fl8*9+T»<  H.  1  if)84934*  T«<  f  .  »*6:>8/  '4,>*  T 
1  •{  0.  Pn.liil  •?.{!'»  TMO.  ft  tiQ.f?41  l  ))))) ) 

RSI i =HS 1 1  *  X 

Y  =  X  /  ?  .  <i 

RSKt  =X»AL00( V  )  »  BS 1 1 *1  .0 

Y  =  Y«  Y 

RSK  i  =HSKl  ♦Y«(H.lS<MlU44*Y»(“n./7?/at>79»YW-t).  1  «156fl9-'  +  Y# 

1  (  -  0.  H  |  9  J  94  0,' ♦  Y»  (-•-  .00  »1  04 H 4*  Y»  ( -n  .  nOfllMhflh  )  )  )  )  )  ) 

RtS!  =RSK1/X 
RQ  TO  AH 


fcCHEB  OHEB 

SUBROUTINE  CHFR ( N I , X 
R I HfcNS ! ON  UX( !  ) 

00  1(1  !=1,N1 
»o  UX(|)  =  0.0 
OX< 1 )  =  1.0 

iiX(?)  =  ? ,  o  *  x 

IF(N1  •  1.  T  .  3  )  RF  TURN 

no  «  i=;i„ni 

ox(  i )  =  ?.n*x«ux<  i-i 

RF  TURN 
FNP 


•  UX) 


>  -UX<! 
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n  r> 


hconsi  roNSi 

BLOCK  DATA 

COMPLEX  A,AA>ANH,CZERO,MASH,AIC  ? 

DIMENSION  AIC(4U,#jn),WASH{4U> 
nOMHON/Cl/Atan). AA(40,8fl>, ANMC 4 0 . 4 0 > , CZERO 

C0MM0N/C?/HKPR(20),ZK6R(?iI)>rHACH(  6 > , FREO< 1 0 ) , N0M«5 > .  1L{50). 

J  HCOR ( 6 ) » ZCOR ( 6 ) *  MXCMN ( 1 1  > ,  WBCN( 1 1 ) . MB I Hi 11 ) , NT (  9 0  >  , 

2  XE<51,YEm,UX<i0>»UYt10>,WXlMM(111,SIX<4U.2>,SCX<l0,2>, 

3  ETA ( 1 1 ) 

COHMQN/C4/CLEN.  NGSKRN, NPY, SOUND, NHACH»HFREG, MAUB, NIOnCX, RHO* 

6  NMODFS#  LCOLL*  LPRWSH, LPRCO, 1 1 V, ! 1 X, NSURF, I SfllAT, FK,FC. 

7  NCOLS,  NOMI Y, MACH ,  XCOLL .YCOLL»PI»U»OWCX#CXMn,I HOD* IRON. 

8  EM.EK»B2»NWlX*NCIX.CRON,NMCY» IFR, El  » E2#  GMY . QWNX, 

9  SN. WRO.NIY,NHCX,NCCX.NWPX#KCPX,NXWING»NYK1nG,NXCS»NYCS 
EQUIVALENCE  (AA.A1C). (A, MASH) 

DATA  Pf /3.1  41  5026n/ 

DATA  HCQR/ 0.0  8566/95  *U.18il  380/9,  ft.  P339569/.0. 23395697, 0.1  8038  ti  79, 

1 fl. 0*566295/ 

DATA  ZOOR/O.n.i  * /6*,  ?a,  (j.  1693953’  ,  0 . 38069U4 1 , 0  •  61  930959 ,  ft  .  83060469  * 

1  0  „9(.6/3476/ 

NOSKRN  SHOULD  RE  COMPATIBLE  WITH  HKER  AND  ZKER  LISTS 
DATA  NfiSKRN/O/ 

DATA  (HKFR<  I  ),  1=1  ,  8  I/O.  05(j6142  7, 0.1 11 19052, 0.1 5685339, 0.181  34  189 
X-ii.  i  81341  89,  n.  15485339,0. 1111905?,  0.05061  497/ 

DATA  <  ZKER ( 1 ), 1=1 ,8 )/U. 019855  07,0. 1  01 66676, 0.2372338  0, 0.4  0828268 
X.  0.591  71  739,  n  .  7  6?  7669  0 , 0 . 89833  <2  4 , 0 . 980 14  49,$/ 

0  •  #**#•****#****#******#<>  ##e***#tf*»*#**o*****#<M*»*«** 

DATA  El  /  0.0 1100001  /,E2/0. 0000  0  «H/.CZERO/(  0.0, 0.0  )/ 

PND 


cxs  xs 

TUITION  XSR.NS.  13,  J3> 

COMM0N/C2/HKER(2n ).ZKER<?0).FHACH(*>.FREQ(1 U),N0H{5) .  IL(5U>, 

1  HCOR  (  6  )  #  ZCOR  ( 6  )  *  WXCHN<  11)/  WBCN ( 1 1  )  . WB I N(11 ) , NT  <  9 ft  > . 

2  XF(l>>,YEt3>»UX(lfl>»UV(lU),MXINN<11),SIX<40.2>,SCX<l'l.2>r 

3  ETAP1) 

COHHON/C4/GLFN.  NGSKRN , NPY »  SOUNH »  NK  ACH,  NF'^EQ,  MAUG,  N!  OvCX ,  RHO, 

6  NHODFS.LCOU.LPRHSH.LPRCO,  1 1 Y#  JiX.NSURF.  !S*L AT , f M,FG, 

/  WOOLS »  NOMI T  #  HACH »  XCOLL • YCQLL, P I »  U.OMCX.CXMw, l  MOD, I  ROM, 

B  EH.EK.82,NWIX,NCIX.CB0N,NHCY, IFR.E1  ,E2,QWY.QWmX, 

9  SN.WBO,N1Y,NWCX,NCCX.NHPX,NCPX,NXWIHO,NYW!nG.NXCS,NYCS 

oo  ro  (in, 4(0, i 

m  GO  TO  (?»,3ll),NS 

20  XS  -  MXC«N<J^)  ♦  W8CN(  J3 )  •  SCXU3.1) 

RETURN 

Jft  XS  -  CXHN  ♦  GOON  *  SCX(I3,?> 

RETURN 

40  GO  TO  (90,6(1). NS 

•>n  XS  =  MX  I HN( J3 )  ♦  HBIN(J3)  «  SIXU3.1) 

RETURN 

XS  =  CXMN  ♦  OB0N  *  SIXU3.2) 

RETURN 

ENP 
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I g^^-zsszmmsive* 


X  I N  F 

FUNCTION  XINT  ( IY# !X.N!Y#NS#WBO»SN) 
r.QHHOH/CZ/U  n  )#XAIC<  10,10,2  )  *  YA  f  CUfl  #  2  >  #B(  40,4ft  )#*<  *0.40)  # 

1  C(4H.4»)#T<4R»4«)#TH<4fl#40)#TR<4fl#4n)#TI(4«#4ii) 

nOHM0N/C?/HKFR(2fl  ) .  ZKER(?»  ) # FHACH(  * ) .  FR6Q<  1  0  ) »  N0H<5  ) ,  it.  <5(1  >» 
f  HCOR( A ) #  ZC0R(6 ) #  MXOHN(  1  i),WBCN(<1  ) ,  WRIK(  1 1  ) #  «T  <  VO  ) , 

i  Xt*(‘>)/Yfe<3>#UX(l0>»UV(1t))»WXIMH(1i>#SIX(40.2)»SCX(10*if)# 

3  ETA  < 1 1 ) 

IF  (NS  .EC.  1 )  00  TO  ?no 
XtNTsWBO#XS<  > »  0 » I  X#  1  ) 

RETURN 

?nn  F F  (YAICdY.i)  .OT.  YE(2)>  GO  TO  ‘3fiO 
XINT=WBO*XS(1 . 1 » IX.l ) 

RETURN 

Tin  SL0PE  =  (WB0»SN»Y(NIY)-YE<2) )/ ( W80«XS(1 » 1  * JX . N! Y >-MBO*vS< 3 #1*  IX# 1 > ) 
XIN1-H0O*XS<1  . 1  .  JX*1  )  +  (YA!CUY.l)-YE<2) ) /SLOPE 
RETURN 
ENR 


*  » 


CARCnOS 

FUNCTION  ARCCGS(X) 

C  •  ••  OFF  I NE  F(X>=A<«  )  +  All >»X«A(?  )»X*»? » . • . »A(ft )*X*»8 
C  THFN  ARCCOS<X)=F<X>Mt-X)**0.5  IF  X.LT.1  AND  .GT.n 
C  *»*  AND  ARCC0S<X)=Pf-(1-A8S(X> )#»«.b*F( ABS(X) )  IF  X.LT.U  AND  .8T.-1 
C  #«*  ARCCOS(X)  IN  RADIANS  *#*  VALID  FOR  0  TO  PI  RADIAN* 

C  »#»  ACCURATE  TO  AT  LEAST  6  SIGNIFICANT  FIGURES  UNLESS  X  APPROACHES  i.tt 
A  0  =  1 .57079631 
Al=-. 21 460184 
A2  =  li.0«904567 
Ai=~. 05072731 
A4  =  (* .  0.1.51  324A 
45  =  -. 02199838 
A6  =  r..  i|1  261238 
A7=-. 00499706 
A8=u .00095128 
IF  (X  .GE.  O.ii)  GO  TO  100 
7=AhS( X ) 

4RCCOS=3. 141  5927- (1  . 0- Z  )  *•« . 5* < AO ♦ A1#Z* A?#Z»»?* A3*Zfr#3 ♦ A4*Z*#4 
l+A6*2*»lJ*A6«Z**64-A7#Z»*7*A8*Z*»ft> 

RETURN 

1  Uli  ARCCOS={1  .n-X)**n.5«(A0*Al«X*A?«X»*2*A3«X««3«A4»X*»4»A5«X**5 
1 *A6*X**6*A7*X*»;*Afl»X*»8) 

RETURN 

END 


< 
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I 


—5  KV**&** 


cmintc 

SUBROUTINE  H!NTC( IS»CNL*CNU#FC) 

0  *##  CHORDw! SE  PRESSURE  INTEGRATION 
n  ###  CNI  =LOHFR  INTEGRATION  LIMIT 
H  •  CNHsUPPER  INTEGRATION  LIMIT 
CL  =  ARCCOS<CNI.  ) 

OU=ARCCOS<CNU> 

SJsHOATC  IS-t  ) 

!F  (IS  .£Q.  ')  00  TO  in 
IF  (IS  .FQ.  o )  00  TO  ?n 

FCS<SIN<<SH  •  .0)*CU ))/(?. (i »(Sl«1.0))~(SlM<(SI-l.0)»Ci»)>/<2.8»CSI 
1  1  •  <1 )  )  - 

2  (S|N((SI*t  .<>)*CL)>/(2.U«(SI*l.ll>)«(SlN((SI*1.0)*Ct  ))/ (?.»•(  Si 

1  i  .  0  )  > 

00  TO  1  (i 0 

in  FC  =  SIN<CU)-CU-SIN<CLHCL 
oo  to  inn 

2D  FC=(SIN(?.0«CU) )/T.«-CU/?.0-(S!N(2.U*CL>)/4.n>CL/2.U 
inn  CONTINUE 
RETURN 
FNB 


i  pji  n  i  5 

r  ...  ?p»!!u?ciNLK,"TSUS'SNt  'swl-'-fsi 

a  ...  iNTMa»nox 

r  **»«*»"«  iwTFOMrioM  S  t 

SL^ARCCOS(SNI  ) 

su*ARrros(SNi/> 

t  is  ,F(J .  r  )  00  ro  Ifl 

US  .go,  ?  J  QO  T0  -,p 

si*fto*rus-i)  G°  T° 

CO  TC  )«»  -'‘•"'(SI»J.o)).(si»((S|..i.o),SU).SIN((^J^'JJ*SLj 

00  TO  1«n  SNU  >**ft.‘>  +  SNL*U.«-SNL*. 

RETURN 


CXLSQ  XLSQ 

SKPROUT ! HE  XLSQ 

HOMPLex  A.AA.AMH.CZERO.KASH/AIC 
HIMFNSION  Air(40,Hn),«ASH(40> 

COHHON/C1 /A(S<i  ),  AAfAO.RQ),  ANM('»('.4rt),CZ&RO 

C0MM0N/C?/HKFR(2n  ),  ZKER(P.)  >.FHACH(M»FR£Q(1  »  )»NGH<5> )  .  IL(5(i ), 

1  HC0RI6 ) » ZC0R(6 ) . WXCMNl 1  1 ).  H8CN(  ll  ) ,  MB  INC  1 1  > ,  W7  ( 90  )  * 

P  Xt(‘>i.YE<-M,UX(lfl>»UY(1ii)*lfXINN(11  >»S!X(4<l.?),SCXUU.fc)« 

•1  F  I A  ( 1 1  J 

COMMON/fM/YI  1  l).XA!C(|tt»1fi.2).YAIC(lO»?)»B(4il»40)«R(4|i,4fl), 
l  C<  1'.  .40  )  .T(4A.4ii)»TH(4  0‘*  * fi > , TR( 40 , 4 u > , T I ( 4  A , 4 0 > 

C0MM0N/C4/ri  FN. NGSKRN.NPY. SOUND .NHACH.N?  RED . MAUG. N ! GwCX , RHO. 

6  NMOnFS.LCnLL.LPRWSH.LPRCO,  IIY,  IIX.NSURF. ISnLAI.FM.FC, 

V  NCOl  S  »  ROM  I  T  ,  MACH. XCQLL  .  YCOLL *  Pi  #  II,  OWCX, CXMw ,  ?  HOD.  IROW, 

H  FM.FK.B %NM!X,NClX>CRON,  NWC  Y » 1FR.E1 . EZ ,  QMY . QMWX, 

y  SN.waO.NIY.NWCX.NCCX.NWPX.HCPX.NXWlNG.NYUIuG.NXCStNYCS 

Foil i valence  <aa,aic>. <a,wash> 
n  =  i 

oo  i =1  *  hCOl  s 

Rl I  =  CABS(AA(  11. m 
If  <RI  I  .LF.EP  )  GO  TO  1.5S 
f  L ( I  )  =  II 

II  r  11  41 

GO  TO  1.16 
1  >s  IL(  I  )  =  “1 

IIP  =  NCOLS  -  1  -  C I - 1 1  ) 

00  » 1  IS  I  1  =  1 1  .  I  I? 

M  l‘>  CA1  I  OORFDl  A  A  (  l  !  «  t .  I  ) ,  I  I ,  {  ♦  1 ) 

MS  noNIFNIlf 

n  SOI  VF  FOR  THF  COFFF I  C  l  FNTS  8Y  BACK  SII8ST  I  IUT  I  ON 

110  II  ■-  WOOLS 

00  •  5 1}  1  =  l.NHOLS 
no  IS.)  L  =  1  »  NMODES 
1'-*»  ANH(I.L)  =  CZFRO 
HO  <  1  -i  J=1  .  NCQI  S 
IF  ML(  1  l  >  .LF.O)  GO  TO  ?1  0 
J!  =  IL(  II  ) 
no  ;  no  i =i,NMonES 
ML  -  WOOLS  «•  L 
irni-wnOLS)  1 7o,  .qi-,p?u 

I/O  IK  11  *  1 

no  i a o  K=1K, NCOLS 

1*m  ANM(II.I)  =  ANM(H.I)  -  A  A  (  J  1  .  K  )  *  ANM(  K  >  L  ) 
ivo  ANM(II.I)  =  (ANHUI. L>  ♦  A  A  (  J  |  .  ML  )  >  /  A  A  (  J  I ,  II) 

Ron  CONI  FNIIF 

Pin  l|  n  -  i 

PPM  PF  TURN 
FNIl 
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CGEOM  GFOM 

SURPOUTINE  GFflM 

COMPLEX  A, AA. ANM.CZFRO.MASH* AIC 
DIMENSION  A1C<4|».  )ft),WASH(4U> 

COHMON'Cl/AlAH  )»AA<4n,ftO)»  ANMf'Wi  .40),CZERO 

COMMOH/C?/HKFR(?n  ),ZKER<?il  >.FMACH(A>.FREO<lO>»NOH(B)  .  I L  ( )» 

1  HC0R(6)»ZC0R(6).WXCNN(l  I  )  *  WBCN(  i  1  ) . WB1NC11 ) , WT(Y0 ) . 

?  XE(‘>),YE(J>,UX(H»>.UY(1«)#WXlHNni>,SlX<4u.,?>.SCX<1U.c>. 

3  ET  A  <  1  n 

COKMON/n.I/Y(i  t  )»  XAICC1  n»1  .<  .2  )  ♦  YAIC(lfl/2)  .BI  4H.4II  )  ,R(  -1U.  40  )  , 

1  C  C  4  ,i .  4  n  ) »  T  ( 4  n ,  4  il ) ,  TH(  4  fl ,  4  0  ) ',  TR  ( 4  0 , 4  ii ) ,  T I  ( 4  i .  4  ii ) 

COMMON/ C  4 /CI.FN .  NGSKRN. NPY »  SOUND. NMACH. NF REQ »  MAUO. N I OmCX  »  PHO » 

(,  NHnilFS.LnOLL.LPRWSN.I  PRCO*  II  Y,  IJX.NSURF.  I  SnL A T . f W, FC  . 

/  NC01.S. NONIT. MACH. XCOLL.  YCOLL  .  P  I  .  II.  QWCX,  CXH*i .  ImoD.  IROW. 

0  EH.FK.B^.NHIX.NCIX.CHON.MHCY. SFR.E1  .F2.0WY  QWwX, 

V  SN- WBO. NI Y.NUCX.NCnX*NWPX»NCPX.NXWlN6.NYU!MG»NXCS*NYrS 

FQII  |  VALFNCE  IAA.AIC).  (A.  WASH) 

C  WBO  =  WING  ROOT  SEMI-CHORD 

C  S  r  semi-span 

n  WTCN  =  WING  TIP  CHORD  -  NORMALIZED  ON  WBO 

C  WTLFN  =  WING  TIP  L.E.  -  NORMAL  I ZFD 

C  SN  =  SFMI-SPAN  -  NORHAI IZED 

T  CBO  =  CONTROL  SEMI-CHORD 

H  FH  :  ?*NWIX+1 

C  FC  =  2*NC I X+ 1 

WBO  =  XF  (vD/^.o 
CLF  N  =  XE ( 4  )  7  WRO 
S  =  YfcCO 

WTTN  =  (  XE  ( -  XE  (  -1 )  )  /  WRO 
WHEN  =  XE  (  2  )  7  WBO 
SN  -  S / WBO 

CBO  =  (  X  F  ( 0  )-XF(4  >  )/?.'! 

FI  =f  W 

F f  -  Fi*P!/V. 

J  =  NW I  X 

COMPOTE  CHORDWISE  INTEGRATION  AND  COLLOCATION  STATION*; 

FIRST  ON  THE  WING  SURFACE 
DO  »>  1  =  1  •  NW  I  X 
F2.  =  F?  -2.«P1 
SIX( J.1  )  =  SINTF2/FJ  ) 

II  t  FLOAT!  J  )/FLOAT(NIONCX)  «  -''.09 
SCXI  I  I  .  1  )  =  -SI X( J. 1  > 

*>  .1  =  J- 1 
r  t  rf  r 

r;  .  f  » •  p  i  /  ? . 

.1  r  NCIX 

r  THFN  ON  THF  CONTROL  SURFACE 

no  ».  I-.  1  .NCIX 
f  ■'  -  I  '*  -■'.•PI 
SI  X  (  J  .  -'  )  =  S  r  N  ( F  9  >  F 1  ) 

II  -  fl  0A  T(|  ) /rLOAT(UIONCX)  ♦  ".9s> 

SCXI  I  i  ,Y)  =  - S  I  X ( J . 2 ) 

J-J-l 

FI  4 • N I  V 

F  1 1 


C  COMPUTE  SPANWISF  INTEGRATION  AND  COLLOCAIION  STATIONS 
DO  i  I = I . NI Y 
III'  =  S  I  N  ( f  *>  /  F 1  ) 

F/  -  f?  *P! 

F T A  (  I  )  =  SIN(F’//Fi  )  112 

H  FV  •-  V9  «P|  ' 


r 


CUHf’tl  I F  WING  SFH  !  -CHORDS  AND  M  l  D  -CHORD  LOCATIONS  AT  THE 
SPANHiSF  COLl OP A  T I  ON  AND  INTEGRATION  STATIONS 
P  I  R  =  YF(  M/YF(.J) 

por  =i.n -Pin 

CRON  =  (XEC>)-YE(  !)}/(?. 0«MBO) 

CXMN  =  CRON  *XF(4)/WR0 
no  •  1=1  »  N  I  Y 

!  F  {  F  T  A  (  I  1.LF.PFB)  GO  TO  1  9 
Fl  '  WTl FN*(FTA( F )-PIH)/POB 
If  (Y(  I  ).!  b.PfH)  GO  TO  1  >1 
F/  -  WTl  f  N*( Y( I  ) —  P I R ) / POS 
GO  10  1 h 


1  ’> 

F  1  -  'l  .  H 

1  s 

1  /  -  H  .  1' 

1  A 

WU  I  N(  I  )  = 

0  .***(}> .  • 

-f  1  ) 

WX I MN(  I  ) 

=  W«1|N(  I  ) 

♦  Fl 

WRCN  f  I  >  = 

n .«»«(?.  ■. 

-r?  > 

1  h 

WXCMN(  I  ) 

=  WRON (  I  ) 

♦  Fv 

R(r TURN 

FN!) 
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iSv.  a. 


CJRAHP 

SUBROUTINE  TR(iMP(NlY»NWCX.NCCX«NXW|NG»NYMlNG»MXCS#NYrS»NIF  . WRU.SN) 
C  *M  TRANSFORMATION  MATRIX  PROGRAM 

r  #«.  TRANSFORMS  Ain  cot  LOCATION  STATIONS  TO  UNSTEADY  AERO  STATIONS 
COHMON'CA/Yt  >  i  >»XAlC(in,1il.2),YAICCin#?)»B(4*  ,4U).R(  Ml . 4>  ) . 

I  C(4..  .  4'M  ,  T  <  M  ,  4.1  > ,  TM(  ,  40  > ,  TR(  40 , 4  |, ) ,  T  1  (  4  >  .  4i,  ) 

i;  **•  7FR0  (TM)  MATRIX  FOR  CHORDWISE  T  RANSFORM  A  T  I  ON 
KROwS=NXWJNG*NYW!NG«NXCS*NYCS 
KCOi S=KROWS 
00  Mi  f =1 , KNOWS 
no  1  tl  Jri  .  KCOLS 
v"  TM{|.J)  =  n.f) 

C  »**  CHORUM  I SF  TRANSFORMATION  (WING) 

IF  (NXW1NG  . FO.  «)  GO  TO  1990 

no  1 fl  M  n  1=1 , NYU  I NG 

CAU  HMAT(NXW1NG»NRSB.NCSB> 

Cam  THAT  (NX  WING.  I .MS1ZE.  I  .  I  . WRO.SN) 

on  0.11  MR=  1  . MSIZF 

no  j  MM  MC=  1  .  NCSB 

TR(MR, MC  )  =« . n 

on  1  DIM  MRC=  '  .MSIZF 

10.M  TR(MR,Mn)=TR(MR.MC)4T(MR,MRC>*M(MRC.MC) 

CAM  CHAT (NWCX . NI  Y. NXW1NG. NYWING. NIF, l . I . NRSC, NCSC. WnO' SN ) 

00  l  II  tl  ^  MR=>  .  HR  SC 
Of)  Ifhi?  MC=  i  .  NCSfl 
T(MR.MC)  =  fl.o 
00  i  n<i?  MRC  -  1  .  NCSC 

1  Cl ■( ?  T<MR.MC)  =  T(MR.MC>  *C{MR.MRC)#TR(MRC.MC) 

KROw=( 1-1 )*NXWINO 

no  i  non  i r= j . nxwing 

LROW  =  KR()Wf  LR 
KCOI  =  ( 1-1 )*NXWf NO 

no  ■  mo  lc=i  .nxwing 
i.coi  =kcoi  » LC 

Idrtn  TM(l  ROW .  LCOL  >  =  T  ( LR .  l.C  ) 
in.if  CON  I  I  NUF 
IRwO  CONTI NUF 

C  *«•  CHORUWISF  TRANSFORMATION  (CONTROL  SURFACF) 

IF  (NXCS  .  FO  •  »  )  GO  TO  P9‘v9 

00  vO.lfl  1=1  .  NYCS 

CAM  HMAT(NXCS.NRSR.NCSB) 

CAM  IMAT(NXCS.  l.MSlZF.i?.  i  .WRO-SN) 

0,i  >'(1(11  MR=  i  . MS l 7F 
DO  :  Mill  MC=  I  .  NCSR 
TR(MR.MC)=0." 

on  >o.ii  mrc-i.msizf 

prnn  ?R(MR.MC)  =  TR(MR, MC ) »  T(MR, MRC  >*M(MRC, MC ) 

CAM  i;MAT(NCCX  .  NI  Y.  NXCS.  NYCS.NH  .P  .  I  .  NRSC.  NCSC  -  W0O.  S'  ) 

00  ;  (III?  MR=  l  .  NRSC 
110  .'(III?  MCsI.NCSH 
T  (  MR.  MC  )  =  H  .  (I 
on  /(III?  MRC  =  1  .NCSC 

Pill*?  T  (  MR  .  MC  )  =  I  (  MR  .  MC  )  »  C (  MR .  MRC  >  * TR (  MRC .  MC  I 
KR(IW  =  NXWIN0»NYWIN(M  (  1-1  )  *  NXC  S 

on  .  linn  lr=  i  -  NXCS 
I  Rnw=  KRdW  *  I  R 

KCOI  'NXWlNO«NYWlN(i*(  F-1  )«NXCS 

on  ponn  I c  =  1 .NXCS 
i  cm  -  Kcoi  ♦Li; 

POOH  TM ( |  ROW . I C0|  )  =  T(LR.| C) 

PHHC  CON M NUF  -114- 


!><)  VO 

r  » *  • 


vn  ,o 

r  •  e  * 


/.(I 

r  » e  * 


U|  .1 


in  f? 


M  / « 


s  n  •<  n 
in  >n 

\  ft  >HI 

\n  i  n 
i  o  g  o 
C  *  » <• 


4  n  ••  i 


-in  r> 


CONTI NUF 

RFARRANfiF  ROWS  ANO  COLUMNS  FOR  SPANWISF  TRANSFORMATION 
C  A I  I  RMAT(NXWING.NYWlNG.NXCS.NYCS.MSIZF) 

do  ;  H‘>n  HR=1  . MS l 7F 
no  ’/»•>«  mc=  i .  ms  I  ze 

T  I  (MR.MDsn.n 

no  ii ‘>«  mrc=-,msjzf 

TI 1MP/MC)=TI (MR,MC)4R(MR.MRC)»TM(MRC,HC) 

7r  pfi  (TM)  MATRIX  FOR  SPANWISF  TRANSFORMATION 
KROwS=N1 Y*<NWCX*NCCX) 

KCOI  S=NXWING»NYWTNG*HXCS*NYCS 
no  (II  1=1  ,  KNOWS 

no  /  ft  .1  =  1 .  kcoi.s 

T  M  (  |  .  J  )  =  o  .  ft 

SPANWISF  TRANSFORMATION  (WING) 

IF  (NYWING  -FO.  0)  GO  TO  4999 

no  'ftnn  I  =1 . nwcx 

CAM  RMATtNYW1NG.NRSB.NCSB) 

CAM  THAT  (NYWING.  ?.MSIZE.  Ml.WRO.SN) 

00  ' ft  *M  MR= 1  . MS  I ZE 
00  .'(mji  MC  =  1  .  NCSB 
TR( MR, MC  )  =  0  .  n 
00  'futl  MRC  =  1  .HSIZE 

TR(MR,MC)=TR(MR,MC)4T(MR»MRC)*R(MRC.MC) 
cam  SMAT(NIY. NYWING. I . NRSS , NCSS . MBO , SN ) 
no  'ftn?  MR= I . NRSS 
no  .inn?  HC  =  l  .  NCSB 
T  (  MR,  MC  )=0.  •• 
no  'ftn?  MRC  = 1  .NCSS 

T(M[/.MC)  =  T(MR.MC>  -  C  (  MR .  MRC  >  *  TR  (  MRC  .  MC  ) 

KROW=(  1-1  )  *  N I Y 
DO  ft  >j n  I.R  =  1  .  N  l  Y 
I.  ROW  =  KROW+LP 
N I T P  =  NX W I NG-  i 
no  .  ft  /ft  J=1  . N I TR 

IF  <  WHO*  X  S  <  i  •  1  .  I  .  i  )  •  L  T .  .*5*  <  XA  IC<  J#  1  ,  1  )  ♦  XA  f  C(  J*1 , 1  ,  M>  >  )GO  TO  'o5u 
CONT inuf 

KCOI  =NYW I NG« (NXWING-1  ) 

GO  10  4n9(l 

KCOI  =NYWING*< J-l » 

no  'ft«n  i  c  =  i . nywing 
i  cm  =  k coi  « ic 

T M ( |  ROW. I COI  )  =  T { IK , 1 C) 

CON  I  I  Nil! 

CONI  I NHF 

SPANWISF  TRANSFORMATION  (CONTROL  SURFACE ) 

II  (NYCS  .  F  U .  *•  )  GO  TO  4999 

no  nun  1=1. NCrx 

CAM  GHAT (NYCS. NRSB.NCSH) 

C  A I  I  TMAT(NYCS.?.MS17E.?.I.WB0.SN) 

(10  .ill  ul  MR=  I  .  MS  I  7F 
no  'in  ul  M  C  =  I  .  N  (’  S  B 
r rtmr. mc  )  =  n .  u 
no  4IMI1  MRC-  ’  .MSI/F 

TR(MR.MO)=TN<MR,MC)«  T(MR.MRC)#B(MRC,MC) 

C  A I  I  SMAKNIY. NYCS.?. NRSS, NCSS.  WRO.SN) 

no  /  n  H  ?  MR= |  NRSS 

no  *11  u?  HC  =  1  .  NCSR 

T  (  MR,  MC  )  =n  .  I’ 

no  -i  II II  ^  MRC  =  '  .  NCSS 

t(mr.mc)  =  t(mr.mc)*C(mr.mrc;^*trcmrc.mc) 


KROW=< 1-1  )*N  J  Y‘ NWCX*  N I Y 
RO  oUMfl  LRsj  ,  N  !  Y 
I.ROW  =  KRO«  +  LR 
NITR=NXCS-1 

no  --n/n  j=i  ,  nitr 

i  f  < W80»xsn  .  i .  i  >  .lt.  jc<  j*i  ,?)*xaic<  j*i  ,1 ,  ■»>  >  >oo  to  «»«s  « 

40  / fl  CONTINUE 

KC0L=NYH1NG*NXMIN(?*NYCS*<NXCS-1  ) 

00  TO  4090 

4  R *>  fi  KCOl  =NYWING*NXW1N0+NYCS*<  J-1  ) 

4 0 9 (i  00  4080  LC=  l  .  NYCS 
LCOI  =  KCOl  HC 

40‘if  TM(!  ROW.LCOL  >  =  T  C I  R,LC> 

4  0  id  CONTINUE 
49V<>  CON  I  I  NUR 

rc  ***  REARRANGE  ROWS  ANll  C0LUHNS  SO  STATIONS  ARE:  STACKED  RnWWISfc 
CAI.  I  RHATTNI  Y.NWCX.NI  Y.NCCX,NS!ZF) 

00  S  0  01  MH= i . NS  I ZE 

no  *» n <n  mc=  '  .  kcols 

TRiMR.MC)=fl." 

no  •  n.n  mrc=  > .nsize 

bn  hi  TR(NR,MC)  =  TR(MR.MC>  +  R(HR»MRC)*TM(MRC,MC) 
no  ‘'049  MR=I  ,  KROWS 
no  MC=1  .  KCOI.S 
Tm ( MR , ho  )  =  n  .  n 
no  MUiR  MRC=1  , KCoLS 

bn  it?  TM( MR »  MC )  =  TM( MR  >  MC )  +  TR( HR#  MRC ) » T I ( HRC»  MC ) 
no  ‘  n^o  j s i , kroms 
no  •  n  >n  j=i , kcols 
TR( I »  J)  =  TM( I .  JT 
b  0  b  n  Tin.  J)sTM(  I.J) 

RETURN 

FNn 
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r  TMA  T  TMaT 

SURROll  T  I  Nfe  I  HA T { NPTS. NO, MS  I ZF, NS.  IY.WBO.SN) 

CQMHON/C  VY(  1  I  ),XA1C(  I  0 , 1 1» .  2  > .  YA  !  C<  1  fi  .  2  >  ,  B  (40 . 4  u  )  ,  R  ( 

1  C<4''.40).T(4n,4!l  ) ,  40  ) ,  TR(  4C  »  40  ) ,  TH4  ‘t ,  4»  ) 

C  ...  GENERATES  (T)*»(-l)  MATRIX 

r.  ...  NPTS  =  NUMBER  OF  AIC  POINTS  ALONG  STRIP  IN  ND  DIRECTION 
c  •  **  MSIZE  =  ORDER  OF  T  MATRIX 

C  *••  NS  =  SURFACE  (1=WING  AND  2=C0NTR0L  SURFACE) 

C  *»*  NO  -  INTERPOLATION  D1RFCTION  U=CHQRDKISE  AND  2=SPAHWI5E> 
IF  ( NPTS  .LT,  4)  MSIZE=NPTS 
IF  (NPTS  .61.  T )  MS  I  ZE  =  3*NPTS-/: 

00  I  J=1 »  MS  I 7E 
00  1  K= 1 »  MS  I ZF 
1  T  (  J .  K  )  =  ii  .  II 

IF  (NPTS  .61.  4)  GO  TO  ^ 0 «| 0 
60  10  (Vf'iifl  if»on,40flt|  ),  NPTS 

r  ...  NPTS=2  (TWO  POINTS  ALONG  STRIP) 
pn  m  T  ( 1  , 1  )  =  1  .  U 
T(?,  1  ) -1  .'i 

IF  (NO  .  Fu.  ')  T(  i  ./>  =  XAIC<1 .  IY.NS) 

IF  (NO  .FQ.  M  T</>»V)  =  XAIC(2»  IY.NS) 

II  (NO  .FQ.  '>)  T { i . 2 )  =  YA ! C ( 1 >  NS ) 

IF  (NO  .FQ.  >)  T(  '.2)=YAIC(2.NS) 

60  10  ft  (i  0  u 

r  ...  *PTS=1  (THRFF  POINTS  ALONG  STRIP) 
i  n  f  T(1  .  1  )  =  1  .0 
T(2.1  )  =  1  .11 
T  <  3 . 1  )  =  1  .  U 

II  (NO  .FO.  >)  60  TO  .■> n  1 0 
C  ...  NPTS=  I  CHORnuiSF  D  I  REFT  l  ON 
T ( 1 , ?  )=XA IC( » . I Y. NS) 

T(i  .  .1  )  =  T  ( 1  ,  2  )  o  *  2 
TC?  .  ?  )  =  X A  ICC> .  t  Y.  NS) 

T  <  ?  ,  3  >  *  T  <  2 , 2  )  *  *  2 
T(T,?)=XAIC(  <. IY.NS) 

T  n .  :n  =  T  (  6  ,  <  )  .  *  2 
60  TO  ft  n  I'll 

r  ...  NPTS=>  SPANWISE  OlRECTION 
(ft  i‘  T(1  .  ?  >  =  Y A  I C (  I  ,  NS) 

T  (  1  ,  T  )  =  T  (  I  »  2  )  »  «  2 
T(?,?)=YAJC( 7 . NS ) 

T  (  *>  .  3  )  r  1  (  2  ,  )  a  a  ? 

TO,2)  =  YAfC(  (.NS) 

T  (.'.')  =  T  (.(.  .  )  «  a  2 
6  0  10  ft  M  ft  II 

r  ...  NPTS-I  (FOUR  POINTS  A! ONO  STRIP) 

<J  M  "  |l  T  (  1  .  1  )  i  I  .  II 
T  (  V  .  1  )  r  1  .  It 
T  (  ( .  1  >  =  I  .  it 
T<4  ,2  )r1  .11 

r  <  *; ,  4 1  r  i .  u 

T  (  ft .  4  )  =  i  .  Ii 
r  (  (  .  4  )  r  -  I  .  II 
r<4  .»>)  =  -).  n 

n  (NO  .FO.  >)  00  TO  4011) 
r  ...  N p T S r 4  FHOROHISF  DIRECTION 
n 1  . 2  )  tXA  IC( I  .  I Y. NS  ) 

T  (i  .  ( i  ?  t  (l ,  ;• )  •  •  2 
T ( 2 , 2  )  t x A  J  C (  ' . IY.NS) 
re*.  ( )  =  T ( 2 , 2  ) * •  2 

T0.2)  =  n.*>MXAIC(  *.  I  Y.NS)^JCU.,|Y»NS>  ) 


bniin 


•>  n  i  o 


T(3,3}s7{  3  »<■•)#*> 

T  O  ■  S  )  =  - 1  (  j ,  r>  j 
T(  >=-T<3.  <  > 

Tf-1 ,  5  )  =  ?»u#|{  \  .  /) 

T(4#rt)r-T(4, 7) 

T  (  %  ,  S  )  s  xA  I C  (  «. JY.NS) 

T(*?.6)rT(5,  »,)<»? 

T{A,r>)sxAIC(4,  I  Y  .  NS  ) 

T(A,6  )  =  7(t>,  >, 

oo  to  #s n n o 

»*  MPTS=1  SPANWISE  direction 
If  T(1 ,?>=YA1C( I , NS ) 

T(1  >3)sTO»y)**i> 

T(? . ? )  =  YA 1C (  1 .  NS ) 

T  (  ?  ,  3  )  =  T  (  ?  >  V* )  *  *  v 

r  (  3 , 3  )  s  T  (  J  ,  y  )  #  #  p 
T(3,  S  )=-T(3,  *> ) 

T  (  3,  ft  )=-7  (3,  , } 

T  ( -i  ,  3  >  =  y  .  |)  *  t  f  < ,  y  ) 

T(4 ,ft)s-T(4,  i) 

T(t),5)  =  YAIC(  ( ,  N  S  ) 

T(fi , b  )  =  yA IC( 1 .NS) 

T(A.M  =  T(6,‘j)*»? 

00  10  A  o  0  u 

*  npts  ,GT.  A 

0  IF  (Nn  *eo*  9>  GO  TO  SSoo 

NPTS  .GT.  4  (CHORnniSF  DIRECTION) 

i(i  .7  )  -  I  »  (j 

HI  .  P  )  =  XA  J  C  (  '  .  I Y , NS ) 

T{1 »3)sT(l>P)#*p 
T<?.1  )  =  ).,! 

T(?,?)=XAIC<?, IY.NS) 

T(HSIZF.HSIZF-*»)s,.1( 

)  =  XAIC‘NPTS,IY.NS) 

r  yp'?  i^,  =  r<MSIZE'MSI7E-J)*^ 

>(MSUF-1  .  MSI  J?F-'M  =  i  .o 

T ( MS  1 /F-1 »  MS  I /F-(  )  =  XAIP(NPTS-i  tv  uci 

NS-Np(s-i'MSr/F,  =  MMS,7K'1 

DO  *  0  I n  N  =  1  , N| 

NRs  ,■  ,  <.N 
NCr  1*N* I 
NP-=N*y 

T  (  NR,  NC  )  =  l.  li 

T ( NR  »  NO  + 1  >  =  XAJC(NP.IY.NS) 

T(NP .NC ♦?)=!( NR »NC*1  ) »  *  p 

N  FrNPTS-3 

DO  ’•«  N=1 , NT 

NRr.  •  N 

NC r  («N-i' 

T(NR.NC)  =  1 
f<NRM  ,NC*1  )r|  .0 
T(NP.NC*3)  =  -i 
T(NR» l ,NC«4 )s>l .o 

T( NR. NC ♦ 4 ) s- T ( NR . NC  +  1  ) 

T<NR.NC*5)=-T(NR.NC«?) 


■  « 


'( /  f> 

r,  **# 

*5  5  II  f| 


S  *5  i  •> 


b  S  /"  ft 

C  *  *  « 

6  0  ii  M 


T  (  NR*  )  .  HC.*?  )  =;  .  ||  •  HNR,  NO  I  ) 

TtNk*} .NR«5)s-T(NR«l .NR»?) 

GO  (0  Mi  Oil 

NPTS  .GT.  4  (SPANWISH  D I RFC  Y ION  > 

T  (  1  .  1  )  =  1  . 1) 

T ( 1  . ?  )  =  YA I C  < 1 > NS  > 

T(1 ,3)sT(t»P)**}> 

T  (  ?  ,  1  )  =  1  .  U 

T(’.?)byA!C(°.MS) 

T(9,  :<)  =  !(?,  £)**?. 

T(MMZF,MSIZF-?)  =  .  .0 
T(HSIZF.MS!ZF-I  >  =  YA I C ( NPTS  *  HS  ) 
T(MSlZF,MSIZF)=T(MSIZE.MS!ZE-l  )«»2 
KMSIZF-1  »  MS  I  7F-?  )  =  i  .1) 

KMSIZF-1  . MSIZF-1  >  =  Y A  I C < NPTS-t . NS > 
KMSIZF-1  .MSI  ZF>  =  T  (MSI  7E-I.HS  I  ZE-1>*»2 

NT=nPTS-a 

00  *  S  i  ll  N  =  1  ,  NT 

NR='  +  <«n 

NC  =  .<«N*  1 

NP  =  N  +  /> 

KNR.NC  )  =  i.i 
T  (  NP  ,  NC*  1  )  =  YA  IT  (  NP.  NS  ) 

T(NR.NC*?)=T(NR.NC*1  )  **2 

NT- NPTS-.1 

no  ►  K  ■> n  N=1 , NT 

NR  =  » »  N 

NC=  *  *  N - 

T(NP,NC)=t  .11 

T(  NRM  ,  NCM  )  =  •  .  0 

T  (  NR  »  NC  ♦.*)  =  -  '  .  n 

T(Nr*  i  .  NCM  ) --1  .n 

KNR.NC*  1  >  =  »  «<YAIC(N*1,NSKYAIC(N*?,NS)) 
KNR.NC*  PKKNR.NC+I  )**2 
T(NR. NC*4 )=-T(NR, NC* 1  ) 

KNR.NC*‘>)  =  -T<NR.NC*?) 

T(NR«  1  .  NC«?)r".:i#  T(NR.NO  1  ) 

T  (  N  R  *  1  »  NC  ♦  S  )  -  -  T  (  NR  *  l  . NC ♦ ?  ) 

INVFRf  T  MAJRIX 
CONTI  NOT 

C  A  |  I.  MINV  (MSI7F.  T  .  C  ) 

RT  TORN 
FNn 
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ccmat 


i  hm 

SUBROUTINE  CMAKNPTS.Nl  Y.NAICPX.NAICPY.NIF  -NS,  I  Y  >  NRS  NCS-WBO.SN) 

C  *<*»  FOR  CHORDWiSF  transformations 

n  *«»  NPTS  =  NUMBER  OF  uHOROWISE  UNSTEADY  AERO  COLIOCATION  STATIONS 
r.  *».  NIY  =  NUMBER  OF  SPANWISf  UNSTEADY  AfcRO  COUOCATION  STATIONS 
r  •  ««  NATCPX  =  NUMBER  OF  CHORDWISE  Alt,  COLLOCATION  STATIONc 
F.  *■>*  NAKPY  =  NUMRFR  OF  SPANWISF  AJC  COLLOCATION  S1AIIONS 
C  IY  r  SPAN  NUHMFR  OF  AIC  STRIP  HtfNfi  TRANSFORMED 

c  NIF  =  COD!  f  OR  |)H  FF  RENT!ATI0N  (|=N0  DERIVATIVE  AND  V-\){  )/DX  ) 

c  **•  MATRIX  SI/F  IS  (NRS. NCS) 

COMMON /C-W  Y(  ’  i  >  .  XA  IC{  I  «> .  I  i.  .  V  )  ,  YA  IC(  ln,P  )  .  KM"  ,  4  ,j  )  ,  R(  i  »  ,  4  11  )  , 

1  C  (  i  >• .  .( >1  )  .  T  t  'I  a  .  4  it  ) »  T  M  <  4  n  .  4  U  )  ,  7  K  (  4  '• .  *  11 )  ,  f  1  ( <•  1 , 4 1.  > 

If  ( NA ICPX  .Of.  .W  00  TO  ^ 

NKS-NA ICPX 
NCSiNA  f CPX 
00  )  |  =  i  , NRS 

00  1  J= 1 . NCS 
1  C  (  I  .  J  )  =  n  .  u 

00  10  inn 
f  NRStNAlCPX 

NCS-  ■'*<NAICPX-'>  ) 

00  n  | = 1 , NRS 
00  ••  J=1  .  NCS 
4  C  (  J  ,,))  =  !'.  11 

1«p  IF  (NCS  .(IT.  M  00  TO  ‘>(in 

IF  (NCS  .EQ.  1.)  00  TO  4fin 

00  10  C  pno  ,  /  .11 ,  ,)•>  I  )  ,  NCS 
C  •  **  Two  POINTS 

?•!»  no  1  1  1  =  1 ,  n 0  irpy 

C  (  I  .  1  )  ~  1  .  ii  J 

C( I.?)=XINT{ 1Y. I.NIY.NS.WBO.SN) 

If  (NIF  .EO.  )  C  (  I  ,  1  )  r  n  .  || 

If  (NIF  .EO.  ')  C(  l,P)=1  .11 

P|C  CONTINUE 
RETURN 

C  »**  TURI F  POINTS 

Tult  DO  1  <i  I  =  I  ,  NA  I  rpx 
C  (  1  , 1  )  =  I  .  (1 

IF  (NIF  .fO.  ')  r  (  I  ,  1  )  r  ij  .  .1 
C(I.?)=XlNT(IY.| , N I Y . NS , WHO . SN ) 

If  (NIF  .  E  0 .  ')  C  (  |  ,  ?  )  =  1  . 11 

C  (  !  .  T  )  =  C  (  |  ,  V  i  0  «  / 

"  -F0>  1  C(  I  .  U=?.i|*XINT(  1Y.  I  .NIY.NS.WHO.SN) 

I  1  »  CON  I  INIK 
Rf  TURN 

C  *«»  four  POINIS 

4‘Hi  DO  1  n  I  =  |  ,  NA  I  CPX 

nx-np  rs- 1 

0(1  '  0/,  Jr  |  ,  NX 

If  (  n  .  ’> »  (  XA  M’  (  J.  I  Y.  NS  )  •  XA  1C(  J4  1  ,  I  y  .  NS  )  )  .0  I  .  X  I  NT(  I  Y.  I  N  I  Y  .  NS »  MdO .  SN 

I  )  )  (.0  TO  -l  II  / 

4  CONI  I  Nil! 

NX-NP  IS 

00  in  4i,h 

4M/  NXr.l 
4"D  KC-. 

If  (NX  .01.  '•  )  K  C  -  4 

C  (  r .  K 1; )  ■  1 .  n 

0 ( I . KC  *  I  >  - X 1 N I ( I Y , I , N I Y , NS  >  WBO . SN ) 

C  (  f  .  K  C  *  >  )  =  C  (  I  .  K  C  *  i  )  *  *  - 
IF  (NIF  .bQ.  I  C<  I  .KC)=07JJ»n 


IF  (  N  |  F  ,  to.  •• )  C<  l  ,  KC«  1  >=1  .  II 

IF  <N|F  .EG.  -)  C(  I .  KC*2  )=v « ‘>»X  I  NT  (  I  Y»  I  »  N I Y , NS »  HBO,  Sm ) 

4 i ft  CONI  I NUF 
RF  TURN 

f  **«  .GT.  FOUR  POINTS 
r>»0  DO  •  1  .1  1  =  1  ,  N  A  I  CPX 
NX=NPTS-1 
00  '  n<>  J=t .  NX 

IF  ( FI .  r'«  (  X  A  1  C  (  J,  l  Y  »  NS )  »  XA  I  C  (  J*  •  *  f-Y .  NS  )  )  •  OT  .  X I  NT<  IY*  1  .  N  I  Y  .  NS .  «B0  *  SN 
1  )  )  <;0  TO  *>fl  / 

CONTI  N II F 

NX  =  NP  rs 
00  10  bu 8 
r>  *i  /  N  X  -  J 

■>‘i»  IF  (NX  .IT.  <  )  GO  TO  ')  II 

IF  (NX  .GT.  NAICPX-;*)  00  TO 
K  (;  -  (  N  X  - ,' )  •  .1  ♦  t 
C  (  I  .  K  C  )  =  1  .  0 

C<I,KC«I)-XINTHY.I  .  NI  Y,  NS.  WHO*  SN) 


C(  1 

,  K  C  ♦  V  )  =  C  (  I 

,KC*. )«*X 

IF 

<  N  I  F  •  F-  0  . 

•>  n<  !  ,KC'  1  >  =  1  .11 

IF 

(  N  f  F  .fcQ. 

i  cd.KC)  =  n.'t 

IF 

(N|F  .FQ. 

.’)  C{|,KC*2)  =  ?.  !*XINT(!Y»  I  » N I Y , NS . WBO . Sm ) 

00 

10  S  1  1- 

C  (  1 

.  1  )  =  1  .  <i 

C(  } 

. P ) s  X 1  NT ( I Y .  I .NIY.NS.WHO.SN) 

C(  I 

.  T  )  r  (;  (  |  ,  t  ) 

e  o  ? 

IF 

(N|F  .FO. 

>  C  (  I  ,  1  )  =  (•.<« 

IF 

( N I F  .FQ. 

’  )  C  (  I  .  ?  )  =  1  . «) 

IF 

<  N  [  F  .  F  Q  . 

•  l  C(  I  ,  1 )  =2.  «»*X  1  NT  (  I  Y  .  1 . NI Y .NS.HHO.SN) 

on 

in  s i  4 

C<  I 

. NCS-P  )  =  1  . 

r.i  i 

. NCS-1  )-> ?NT( I Y.  I  ,  N! Y,  NS, WHO,  SN) 

C(  I 

. NCS )-C<  1  . 

NCS-  !)••/' 

IF 

<N|F  .  E  0 . 

'  >  n <  i  .NCS-? )  =  «.(» 

1  F 

(N|F  .FO. 

1  C(  1  .  NCS  -  1  )  =  \  .It 

IF 

( N l F  . FQ. 

'  )  c  (  I  ,NnS)  =  /.ll*XlNT<  IY.  I.NIY.NS.WBO.SN) 

!>  ,  <•  CON  i  I  Nllf 
RFTllRN 

F'NO 
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r  ';HfcT  jM/.  5 

SiiflPOlJT  1  SM  A  T  (NIY.NAIGPY  •  NS »  NRS  »NCS#WHO*SN) 
r:  •  SPANWISf  IRANSfORMAT  ION 

r,  *»«  M  |  Y  =  NlJMHFR  or  SPANWISF  UNSTEADY  AtRO  COL  I  00  A  T I  ON  STATIONS 
r  NA1<  PY  =  NIIMRFR  Of  SPANWISF  A 1  •*  001.  LOG  A  T  JON  STATIONS 

C  **»  NS  -  SURFaCI  <1=WING  AND  P=CONTROL  SlfRf  ACF  ) 
r  #«*  MATRIX  SIZF  ?  S  NRS  BY  NCS 

COMMON  /(’.</  Y  <  1  >  )  .  Jt  A  I  C(  I  A  » t  i.if>»YAlC(l(*»?)*B(4,l»4u)»R(  *')••<  ’  ) » 

1  c<  1  n  .  4n  )  .  T  I  in  ,  4.1  )  ,  TH(  4  0  ,  4  0  )  ,  TR(  4I> ,  4|»  )  ,  T  I  (  A"  ,  4n  ) 

IF  (NAICPY  -OT.  3)  GO  TO  M 
NRS - N I Y 
ncs=na ICPY 
00  /  1  =  1.  NRS 

00  c  J= I • NCS 
'  C(?.J)  =  1.  .!> 

CO  10  ’  ffi 
••  NPS--NIY 

NCS-  1 « ( N A ICPV-  1  ) 

00  ^  1  =  1*  NRS 
00  *.  .1  =  1  •  NCS 
0  C  <  !  .  J  )  =  I. .  i 

1  IF  ( NCS  .  0  T .  .)  fiO  10  Sim 
IF  (NCS  .FQ.  1  00  TO  'Min 
00  10  (  Pint ,  /'i.i  .  .Vi  •.  > ,  NCS 
C  ***  TWO  POINTS 
?<  n  no  /fin  l  =  1 .  n  i  y 

C  (  I  ,  1  )  =  1  . 1) 

CM.? )  =  WRU*SN*Y(  I  ) 

P'."  CON  I  1  Nllf 
RETURN 

C  *»*  THPI  F  POINTS 
Suf  00  fi  I  I  =  1 . N I Y 
C  (  I  .  1  )  =  I  .  i. 

C<  I  .  P  )  =  WR (1* S N  *  Y  (  5  > 

c  ( i .  =  c  ( I . :  i » • ' 

U>'  CONTINUE 
RFTI'RN 

C  **'  FOUR  POINTS 
4.  (i  0  0  '.Oil  1  =  1  ,N?  Y 
M,-' 

1 1  (  WMO#SN»  Y  (  I  1  .IT.  >».».*(  YAICIP  .  NS)«  YA1C(  1,NS  )>  >  ICr) 

'*(1.10=1.0 
” (  I  .  I  C ♦ 1  )  =WPO#SN» Y ( I > 

c  <  i .  i  c  •  ? )  =  c  <  r .  i  c  *  i  > » •  p 

4  CONI  I  Nllf 
PI-TURN 

r  •  ,<;T,  FOtlR  POINTS 

*>  •  '*  110  •  '>  i  i  =  J  .  N  I  Y 

N |  -NA |CPY-P 
0  0  •.?•»  .)  =  I  .  N  I 

U  ( It  .•>»  (  YA  |  P(  I,  NS  )  ♦  YA  IC(  JM  .  NS  )  )  .OT.  WRU*SN«  Y  (  I  )  )  nO  10  02M 
'»/•>  CON  f  I  NIK 

rr;-/*NA  icpy-'I 
0(1  (0  *./4 

'  1  I  c-  (  J-P  )#.)♦  i 

If  (  .1  .If.  11  IC  =  I 
/  4  o  i  i .  i  o )  =  i .  o 

Cl  f  .  1 1>  )  )=WDO  =  SN*Y(  I  ) 
c  (  I  ,  I  C  *  /  >  =  C  (  I  .  I  C  ►  •  )  *  *  • 
b  '■  CON  I  I  Nllf 
MR  TURN 
F  NP 
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C  °N  4  T 

C  »  »  o 
r  »  a  » 
C  S  0  o 


1  'III 


0  I, 


1  .  •> 
1  •' 


RHaT 

SUOPOIITINJ-  RM  A  T  (  NXW I  MR,  NYW I  NO  •  NXCS,  NYCS .  MS  !  ZE  ) 

REARRANGES  Ain  STATIONS  INTO  PROPER  SFQOFNCF  I  OR  SPAmWISF 
INTI  RPOI  ATION 

MATRIX  SIZE  IS  MSlZF=NXWf NG«NYWINfi»NXCS*NYC5  { SQUARE  ) 

COMMON 'C-W  Y(  *  .  ),XA!C<!f>,t'l.2),YAlC(in,?),B(4l',4|j),R<  ui.AO), 
CD  ..4<i).T(/jn>4ii)<TM(4|'*»4u),TR</in,A(i),TI(4  • ,  4h  ) 

C?FPO  =  ll .  I' 

MSI  /  T  =  NXH|Nf:«NYMINR*NXrS*NYCS 
no  -(Ml  J  =  l  .  MS  I  7E 
no  •  I'  l  .j-  1  ,  MSI  7E 
R<I.J)=r7FRn 

II  (NXWINO  .M).  ")  OO  TO  /">» 

K  -  1 
KK" 

I  l =  NYU  I NO«NXW I  NO 
no  n  .i  1  =  1,11 
R  (  I  .  K  )  =  l  .  'I 
K=K«NXHING 

If  (K  . 0  T .  II)  KK  =KK  *  I 
If  (K  .OT.  ID  K  =  KK 
CONllNOf 
CON  I  I  NOP 

If  (NXCS  . F 0 .  i)  00  TO  .**>•! 

I  I =NXCS«NYCS 
K=NyWING*NYWING» I 
KK -t.XU  I  NR»NYW  I  NO  ♦  t 
no  n |  =  |,II 

I  Kr  |  ♦ NXW I  NO  »  NYU  I  NO 
R  (  I  Y  ,  K  )  =  )  .  I' 

K  -  K  « NXCS 

If  rK  .01.  MS  I 7F )  K  K  =  K  K ♦ 1 

II  (K  .OT.  MSI  7  F  )  K  =  X  K 
CON  I  I  NO! 

CON]  I  NIJfc 
Rf  TORN 
FNI) 
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fHMAT  RMaT 

StIKpnilT  I  Nl-  DMA  T( NPTS ■  IRS'  I  l-S  ) 

coupon m  w  y  ( 1 1  > ,  xaii;(  t  n ,  i  >1 .2 ) .  ya  icn  p,?  >  .bmii  »4ij  )  ,R<  •»  •  ) . 

J  C(  * .  1  »  )  .  T  {  ft ,  4  n  .  ^  II  ) ,  TRMft  ,  -ill  )  ,  T  I  (  a  • ,  Ail  ) 

f*  R  -  P,{  |«S*  ICS)  HA  JR  IX 

C  NPTS  =  NUMBFR  nr  AlC  STATION  ALONG  STRIP  (FITHFK  CHOnURISfc  or 

r  »**  SPANHISF) 

ICS=NPTS 

IF  (NPTS  .  G  1  .  « )  GO  TO  /On 

I  RS^NPTS 

no  •  n  i=i , ips 
no  *-n  j=i ,  ms 

R  (  I  .  J  )  =  O  .  II 

IF  (I  .  FO.  J)  R  <  I  •  J  )  =  i  .  II 
>»  nONUNIIF 
Rf  TURN 

pill)  1RS  =  <S.  (NPTS-i  )«  > 
no  mi  1=1 , IRS 
no  ‘ii it  j=i  ,  ms 
J  id  R  <  I  .  J  )  =  H  .  H 
im  .i  >  =  i .  ii 
r  m .  p  >  =  i .  >■ 
r<  irs.  h:S)  =  i 
R ( IPS-I . ICS- '  >  *  l  .  * 

If  ( NPTS  .EO.  * )  GO  TO 
K=NPTS-4 

RO  'H.i  1  =  1,  K 
NR  =  > •  \*  1 
NHr:  ♦  j 

Tjii  R(NP.NC)  =  1  .n 
4  in  Rl TORN 
FND 


I 


1 


-124- 


FMINV  H I N V 

SUP f-'OllT  1  Nh  HINy  (NM.A.1I) 

0  |  MF-NS  !  ON  A(4<i.4»i),U(4I',4ii) 

0(1  •  <>:i1  1  =  1.  NH 

00  Mi.n  J  =  1  .  NM 
II (  I  .  J)s«.i» 

U  (  I  .FO.  J)  tl<  !  »  J  )  =  I  .  » 

Vftnl  C0N71NUF 

F  PS  id.  il  MM  Hun  i 
no  '■iih  l  =  i  ,nm 

K  =  I 

1 1  (I  -NM  )  Oi.*7  .  .  9H  i  /  ,  9ii?1 
mi  '  i  n  <A(i,i)-fPS)  y..ii‘.  .-/'mi  a.  on  it/ 

Ml  '.  II  ( -A  (  I  »i  )-H»q  )  <  ft  u  ft  ,  0  |j  ft  f»  ,  9  n  II  7 

Ml'ift  K  :  k  ♦  i 

oo  mum  j  =  i.nm 

U(  I  .  J )  =  U c  I  »  J  )  ♦  II  ( K  .  J  ) 

9ft A( I , J)=A( I , J)'A(K.J) 

00  10  9  H  V  ! 

90  .7  01 V  =  A(  1,1) 

no  ft  II  li «  J  =  1,NM 

|}<  l .  J)  =  II<  I ,  jwniv 
yni o  a ( i . j )  =  a ( l , ji/oiv 

00  Ml  '  MM=  i  .  NM 
Oft  T=A(MM,  I  ) 

IF  (AMStPELD-FPS)  vfl1‘-.9il1*.9«1ft 

ftnjA  ff  (HM-I)  9-»  * -l  ,9ft  .  ^,9.|1  « 

0-1  I  n  00  ft 0  I  1  J  =  1  ,  NM 

IJ(  MM ,  J )  =  0 ( MM .  J ) -  U(  I  #  J ) *DEL T 
9  0  i  i  A  (MM.  J)=A  (MM.  .))-A(  I  ,  J)*0F!.T 
90(9  CONTI  NIIF 

1)0  <vft  M  I  =1  ,  NM 
00  '<0  M  J  =  1  ,  NM 
9  0  i  <  A  (  I  .  J  )  =  U I  I  .  J  > 

RETURN 


FART  IV 


SECTION  B5.0 


FIDW  CHARTS  FOR  SUBSONIC 
AIC  COMPUTER  PROGRAM 


I 


dvmr*  v  m. 

I  5  I  ( — 


JTOWTHR  To  SIWKTIVEL- 

J  1 

“nj  VRn^{hfhfi) 

\  J 

TO**\r2.(bP|/FW  j 

c*r*-2.o«pi/FC 


nKPf- rr  Tn  soo 

R^R 


MVOUMHM 

»>pfhio«m¥.w 

l-sB*»UWM> 


UnVW?  Tn  M1JO(T|VK 
K*T 
i 


,  ,  TO- W>H  TO  St  lltX  Tl 

»««i..TK.o  LacsJ  kxt 

2 


ftKPKTf  m  mo 

R* 

fffi-I.IM . ■*' 


»)  i  [-US-J  ninj(  ihikmihtmimom /«»nt»2.o*pi  > 


Mwmc 

V1«KS-*NMMMfcV*  »l  >«•  VltTOW  (S 


Dl  HK  V  J  ICNEO  VAWI4H.I.  KS 


<■>  MW.  simvrc 


«V«W.  STW-VW 

<«np  2.3 


SV'tW.  PTOAMCS 
iro  i 


STORM*? 
C  SO 


SV*W.  «TORV»S 


( 
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.MIKXTlMv  (TXT  Cl  TO) 


2MK  • ' T  M I 20  H< ,  V -  HAKE,  26H/ 
orr  xcfiRfc  it  i)  ,i-i  ,3)/*h*(  sc,  s*htai  hh*i  v:*tai  \j 


HV|'lTH5  «1  TO 
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